8  z  s  i  z  r  v  ov  A3oo  m  siic 


i 


NORDA  Report  22 


<  TEST  AND  EVALUATION  OF 

CGC  POLAR  STAR 
WAGB  10 

3 

VOLUME  III  BACKGROUND 

JAMES  P.  WELSH,  JR.,  EDITOR 


Polar  Oceanography  Branch 
Oceanography  Division 
Naval  Oceanographic  Laboratory 

September  1978 


Prepared  for: 

Coast  Guard  Research 
and  Development  Center, 
Groton,  CT. 


Approved  for 
Public  Relee— 


vr- 

V  • 
p?  \n 


W  1  6  1982 

A 


NAVAL  OCEAN  RESEARCH  AND  DEVELOPMENT  ACTIVITY 

NSTL  Station,  Mississippi  39529 


82  11  15  071 


L 


TEST  AND  EVALUATION  OF 
CGC  POLAR  STAR 
WAGB  10 


VOLUME  III  BACKGROUND 

JAMES  P.  WELSH,  JR.,  EDITOR 

Polar  Oceanography  Branch 
Oceanography  Division 
Naval  Oceanographic  Laboratory 

September  1978 


Prepared  for: 

Coast  Guard  Research 
and  Development  Center, 
Groton,  CT. 


Approved  for 
Public  Release 


NAVAL  OCEAN  RESEARCH  AND  DEVELOPMENT  ACTIVITY 
NSTL  Station,  Mississippi  39529 


FOREWORD 


f 

CGC  POLAR  STAR  is  the  first  of  a  new  class  of  American  icebreakers  built  for 
and  operated  by  the  U.S.  Coast  Guard.  The  vessel's  design  incorporates  many  new 
and  sophisticated  systems  not  found  on  earlier  icebreakers. 

The  primary  objective  of  the  test  and  evaluation  program  was  to  examine  the 
performance  characteristics  under  actual  polar  operations.  Two  field  trials  were 
conducted.  The  first  in  the  Arctic  (Appendix  A)  and  the  second  in  the  Antarctic. 
During  these  field  trials  approximately  120  individual  parameters  were  measured  and 
recorded  on  magnetic  tape  (Appendix  B).  Additionally,  the  physical  properties  of 
the  sea  ice  in  which  the  ship  was  operating  were  measured. 

I  Documentation  of  the  complete  program  including  preliminary  screening  of  the 

.  data  has  been  accomplished  by  NORDA  under  contract  to  the  Coast  Guard  Research  and 

Development  Center.  The  documentation  consists  of  the  following  four  volumes:  I. 

\  Antarctic  trials,  II.  Test  Plans,  III.  Background,  and  IV.  Instrumentation  Manual. 


•  1 


EXECUTIVE  SUMMARY 


H 


Section  I  is  a  description  of  the  POLAR  class  vessels.  The  description  was  pre¬ 
pared  by  excerpting  the  work  of  various  individuals,  including  CWO  V.B.  Midgette  and 
LCDR  T.E.  Braithwaite  and  the  Welcome  Aboard  Pamphlet.  This  information  has  been 
published  elsewhere  in  Coast  Guard  documents  but  is  reproduced  here  for  ready  ref¬ 
erence. 

Section  II  through  VI  are  the  work  of  CDR  G.P.  Vance,  USCG.  During  the  summer 
of  1975  CDR  Vance  attempted  to  predict  the  icebreaking  performance  of  the  POLAR 
class  vessels  prior  to  their  first  use  in  the  ice.  This  was  done  to  assist  in  the 
development  of  the  icebreaking  test  plan  for  these  vessels.  He  based  his  predictions 
on  the  results  of  the  various  equations  and  theories  available  in  the  literature  at 
that  time.  Vance  evaluated  the  adequacy  of  these  equations  and  theories  by  comparing 
the  results  of  other  classes  of  icebreakers  to  the  actual  performance  of  those 
vessels.  The  magnitude  of  measurement  errors  for  various  parameters  and  physical 
properties  during  the  testing  were  estimated  to  determine  their  impact  on  the  pre¬ 
cision  of  the  overall  test  results. 


RALPH  R.  GOODMAN 
TECHNICAL  DIRECTOR 
NORDA 


i 


CONTENTS 


Page 


LIST  OF  ILLUSTRATIONS  vi 

LIST  OF  TABLES  vi 

I.  GENERAL  COMMENTS  ON  POLAR  STAR  1 

II.  COMMENTS  ON  VARIOUS  METHODS  OF  EVALUATING  AND/OR  PREDICTING 

ICEBREAKER  PERFORMANCE  BY  CDR  G.  P.  VANCE  2 

A.  Introduction  2 

B.  Category  I  3 

1.  Discussion  3 

2.  Slush  5 

3.  Ramming  5 

4.  Extraction  5 

C.  Category  II  5 

D.  Summary  6 

III.  PREDICTED  OPEN  WATER  PERFORMANCE  CHARACTERISTICS  OF  THE 

POLAR  STAR  BY  CDR  G.  P.  VANCE  6 

A.  Abstract  6 

B.  Basic  Equations  6 

C.  Expected  Values  of  Thrust  11 

D.  Expected  Values  of  Shaft  Torque  11 

E.  Example  of  Torque  Calculation  12 

F.  Propeller  Performance  12 

IV.  PREDICTED  ICEBREAKING  PERFORMANCE  CHARACTERISTICS  OF  THE 

POLAR  STAR  BY  CDR  G.  P.  VANCE  12 

A.  Abstract  12 

B.  Introduction  27 

C.  Nomenclature  for  Ice  Resistance  Program  (Input)  27 

D.  Nomenclature  for  Ice  Resistance  Program  (Program  Constants)  31 

E.  Nomenclature  for  Ice  Resistance  Program  (Output)  32 

F.  POLAR  STAR  Station  Beam  Coefficient  32 

G.  POLAR  STAR  Station  Area  Coefficient  33 


iii 


9 


1 


Page 


J 


N 


V. 


i 


« 


H. 

Transverse  Spread  Angle  Complement  [GBETA  (I)] 

34 

I. 

Station  Waterline  Inclination  [GDELT  (I)] 

36 

J. 

Expected  Thrust  Values  when  Breaking  Ice 

36 

K. 

Expected  Torque  Values  when  Breaking  Ice 

37 

L. 

Calculation  of  Expected  Ice  Resistance  Torque 

37 

M. 

Prediction  Program  Verification 

38 

PREDICTION  OF  ERROR  SOURCES  FOR  FULL  SCALE  TEST  OF 

POLAR  STAR  BY  G.  P.  VANCE 

38 

A. 

Abstract 

38 

B. 

Introduction 

41 

C. 

Program  Parameters 

41 

1.  Milano's  Algorithm 

41 

2.  Vance's  Algorithm 

41 

D. 

Effect  of  Thickness 

42 

1.  Run  21 

42 

2.  Run  30 

43 

E. 

Summary  of  Thickness 

45 

F. 

Effect  of  Flexural  Strength 

45 

1.  Run  32 

45 

2.  Run  21 

45 

3.  Run  30 

45 

4.  Run  27 

45 

5.  Run  36 

46 

6.  Run  37 

46 

G. 

Summary  of  Flexural  Strength 

46 

H. 

Effect  of  Ice  Specific  Gravity 

47 

1.  Run  28 

47 

2.  Run  21 

47 

3.  Run  33 

47 

I. 

Summary  of  Specific  Gravity 

47 

J. 

Effect  of  Velocity 

1.  Run  21 

48 

2.  Run  30 

48 

« 


iv 


K.  Summary  of  Velocity  Effect 

L.  Effect  of  Compressive  Strength 

1.  Run  35 

2.  Run  21 

3.  Run  34 

4.  Run  17 

M.  Summary  of  Compressive  Strength  Effect 

N.  Effect  of  Friction  (4) 

1.  Run  39 

2.  Run  21 

3.  Run  40 

4.  Run  41 

O.  Summary  of  Friction  Effect 

P.  Effect  of  Propulsion  Equipment  Measurements 

Q.  Equations 

1.  Effect  of  Thrust 

2.  Effect  of  Torque 

R.  Electrical  Power  Measurement 

S.  Summary  of  Propulsion  Equipment 

T.  Effect  of  Draft 

1.  Run  45 

2.  Run  44 

U.  Summary  of  Draft  Effect 

RAMMING  AND  EXTRACTION  PREDICTIONS  FOR  THE  POLAR  STAR 
BY  G.  P.  VANCE 

A.  Abstract 

B.  Introduction 

C.  Ramming  Capability 

D.  Extraction  Difficulty 
REFERENCES 


Page 


APPENDICES 


Appendix  A: 

Appendix  B: 
Appendix  C: 
Appendix  D: 


Curves  of  Resistance  versus  Velocity  and 
Thrust  versus  Velocity 

POLAR  STAR  30  ft.  Draft  Particulars 

Sample  Calculations 

Literature  Citations  and  Selected  Bibliography 


63 

71 

73 

81 


vi 


LIST  OF  ILLUSTRATIONS 


Page 


Figure  1: 
Figure  2: 
Figure  3: 
Figure  4: 
Figure  5: 
Figure  6: 
Figure  7: 
Figure  8: 
Figure  9: 

Figure  10 
Figure  11 
Figure  12 
Figure  13 
Figure  14 
Figure  15 
Figure  16 
Figure  17 
Figure  18 


POLAR  STAR  Open  Water  Resistance 

Open  Water  Thrust  Torque  Curves  (POLAR  STAR) 

Stock  Propeller  Performance  on  POLAR  STAR 

POLAR  STAR  Tow  Rope  Thrust 

POLAR  STAR  Section  Area  Curve  (28  ft.  W.  L.) 

POLAR  STAR  Run  21 
MACKINAW  Run  24 

POLAR  STAR  Icebreaking  28-30  ft.  Drafts 
Prediction  of  Thickness  of  Standard  Sea  Ice  Broken  as  a 
Result  of  Ramming  as  a  Function  of  Impact  Velocity 
Extraction  Difficulty  versus  Impact  Velocity 
POLAR  STAR  Run  2 
POLAR  STAR  Run  6 
POLAR  STAR  Run  7 
POLAR  STAR  Run  10 
POLAR  STAR  Run  14 
POLAR  STAR  Run  15 
POLAR  STAR  Run  28 
POLAR  STAR  Run  30 


LIST  OF  TABLES 


Table  1: 
Table  2: 
Table  3: 

Table  4: 


Table  5: 
Table  6: 
Table  7: 

Table  8: 
Table  9: 
Table  10: 
Table  11: 


Summary  of  I/B  Prediction  Techniques 
POLAR  STAR  Open  Water  Performance  (Test  6-A) 

The  Test  Program  for  the  U.  S.  Coast  Guard  Proposed 
Icebreaker  Design  (M-14-3) 

A  Brief  Description  of  the  Stock  Propellers  Used  in  the 
Testing  Program  for  Simulating  Ahead  Modes  of  Operation 
for  the  U.  S.  Coast  Guard  Proposed  Icebreaker  Design 
(M-14-3) 

Predicted  Free-route  Propulsion  Coefficients 
Predicted  Tow-rope  Pull  Propulsion  Coefficients 
Summary  of  Predicted  Towing  Performance  for  the  M-14-3 
Design  Icebreaker  Shaft  Horsepower-20,000  SHP/SHAFT 
Summary  of  Runs 

POLAR  STAR  Input  Data  ( ICEDTA)  28  ft.  Draft 
WINDCLASS  Input  Data  (WINDTA) 

MACKINAW  Input  Data  (MACDTA) 


I.  GENERAL  COMMENTS  ON  POLAR  STAR 

The  USCGC  POLAR  STAR  (WAGB  10)  is  the  first  new  icebreaker  to  be  built  for 
the  United  States  since  1954.  The  Lockheed  Shipbuilding  and  Construction  Company 
of  Seattle,  Washington,  built  the  POLAR  STAR  and  is  presently  completing  her  sister 
ship,  USCGC  POLAR  SEA  (WAGB  11).  Eventually  these  two  ships  will  operate  out  of  a 
support  complex  located  in  their  home  port  of  Seattle,  Washington. 

POLAR  STAR'S  design  resulted  from  three  years  of  research  and  testing.  The 
design  incorporates  a  number  of  improved  and  innovative  features  that  affect  nearly 
every  aspect  of  her  operation.  Much  of  the  equipment  installed  aboard  is  highly 
sophisticated  such  as  the  propulsion  control  system  and  electronic  navigation  equip 
ment.  There  are  three  computers,  one  for  monitoring  engineering  machinery,  one  for 
handling  oceanographic  data,  and  one  for  navigational  capabilities.  In  keeping 
with  current  technology  and  the  desire  to  reduce  manpower  requirements,  automation 
fulfills  functions  to  a  degree  unprecedented  in  the  Coast  Guard. 

The  hull  shape  for  POLAR  STAR  is  derived  partially  from  studies  made  by  Capt. 
Roderick  M.  White,  USCG,  whose  ideas  were  incorporated  into  the  bow  design  of  the 
icebreaking  tanker  MANHATTAN.  Basically,  the  hull  is  configured  to  maximize  ice¬ 
breaking  capability  through  the  combined  forces  of  the  ship's  forward  motion  (as  it 
rides  up  onto  the  ice)  and  the  downward  pull  of  gravity. 

The  shell  plating  and  support  structures  are  fabricated  from  A-537  steel  hav¬ 
ing  chemical  and  mechanical  properties  which  can  withstand  heavy  impacts  and  load¬ 
ing  at  low  temperatures.  Essentially  three  different  types  of  steel  have  been  used 
in  the  construction  of  the  ship  (A-537,  HY-80  and  A- 131).  The  icebelt  plating  in 
the  forward  and  after  portions  of  the  hull  is  1  3/4"  thick.  This  shell  plating 
thins  down  to  l!a"  amidships.  Portions  of  the  hull  that  are  above  the  waterline  and 
clear  from  ice  are  V  thick. 

The  propulsion  system  centers  on  a  diesel-electric/gas  turbine  plant  that  pow 
ers  three  shafts.  The  diesel  generators  supply  AC  electrical  power  which  is  recti¬ 
fied  for  the  DC  propulsion  motors.  The  gas  turbines  are  coupled  via  a  reduction 
gear  directly  to  the  shaft.  Either  the  diesel-electric  or  the  gas  turbine  mode  of 
operation  can  be  used,  but  not  both  together.  The  diesel -electric  units  produce  a 
total  of  18,000  continuous  shaft  horsepower  and  are  the  normal  mode  of  propulsion. 
The  gas  turbine  units  produce  a  total  of  60,000  continuous  shaft  horsepower  for  es¬ 
pecially  heavy  icebreaking  or  emergencies.  Due  to  the  turbines'  higher  rate  of  fu¬ 
el  consumption,  they  will  be  used  sparingly. 

The  propulsion  machinery  is  spread  over  five  separate  compartments.  The  for¬ 
ward  two  compartments  house  diesel  generators,  supplying  electricity  for  propulsion 
and  shipboard  use.  The  third  engineroom  aft  contains  the  gas  turbines.  The  after 
two  compartments  contain  reduction  gears  and  electric  motors. 

The  three  shafts  (34"  diameter)  turn  four-bladed,  controllable  pitch  propel¬ 
lers  having  16'  diameters.  This  type  of  propeller  permits  almost  instaneous  con¬ 
trol  of  the  ship's  direction  of  propulsion.  Increased  maneuverability  is  possible 
with  this  multiple  shaft  configuration.  Additionally,  the  single  rudder  being 
mounted  directly  aft  of  the  center  line  propeller  will  afford  greater  maneuverabil¬ 
ity  than  that  of  other  U.S.  icebreakers. 

The  shafts  and  motors  are  non-reversing  and  maneuvering  control  is  obtained 
by  adjusting  the  pitch  of  the  blades  on  each  propeller.  The  continuously  turning 
shafts  should  be  valuable  during  icebreaking  as  studies  indicate  an  icebreaker's 
propeller  it  nost  vu1'  .'able  when  stopped  for  reversing  direction.  With  a  control¬ 
lable  pitch  ^.^’ll'  the  need  for  reversing  the  revolution  of  the  shaft  is  elimi¬ 
nated. 


A  computer  controlled  monitoring  system  keeps  watch  over  the  equipment  in  the 
engineroom.  Selected  machinery  parameters  are  continuously  scanned  and  recorded. 

If  critical  values  are  reached,  the  monitor  system  sounds  an  alarm  and  automatically 
secures  certain  equipment  before  permanent  damage  is  done. 

It  is  impossible  to  describe  even  briefly  all  of  the  equipment  aboard  POLAR 
STAR.  However,  the  following  items  may  be  of  special  interest. 

One  innovation,  as  far  as  the  Coast  Guard  is  concerned,  is  the  central  hydrau¬ 
lic  system  which  provides  a  single  hydraulic  power  source  for  the  entire  ship.  This 
avoids  the  necessity  of  generating  the  hydraulic  power  electrically  at  each  individ¬ 
ual  machinery  location.  The  15  ton  cranes  aft  and  3  ton  crane  forward,  for  example, 
are  both  operated  off  this  central  system. 

There  is  a  heeling  system  aboard  utilized  to  rock  the  ship  as  an  aid  in  free¬ 
ing  herself  when  stuck  in  ice.  The  system  consists  of  three  pairs  of  tanks  connect¬ 
ed  by  flumes  with  axial  flow  reversible  pumps.  Operated  hydraulically,  the  pumps 
transfer  the  contents  of  each  (approximately  35,000  gallons)  or  all  three  tanks  to 
its  mate  on  the  opposite  side  of  the  ship.  This  is  accomplished  in  50  seconds  gen¬ 
erating  24,000  foot  tons  of  torque  on  the  ship. 

Icebreakers  are  infamous  for  their  heavy  roll  in  the  open  sea.  To  dampen 
this  effect  a  passive  roll  stabilizing  system  is  installed.  A  continuous  port  to 
starboard  tank  is  kept  approximately  40a  full  of  water.  The  placement  of  vertical 
pipes  within  the  tank  produces  a  nozzle  effect  to  keep  the  "sloshing"  water  out  of 
phase  with  the  roll  of  the  ship,  thus  tending  to  counteract  the  roll.  The  effect 
is  a  dynamic  one,  that  is,  it  requires  the  ship  to  be  rolling  before  it  operates. 

One  special  feature  of  POLAR  STAR  is  the  aloft  conning  station  approximately 
104'  above  the  waterline.  This  platform  provides  a  necessary  valuable  visual  exten¬ 
sion  to  the  conning  officer  as  he  carefully  searches  for  the  best  route  through  con¬ 
centrated  ice  fields.  Complete  rudder  and  engine  controls  and  necessary  navigation¬ 
al  devices  are  located  in  the  station  aloft.  Access  is  via  a  ladder  inside  the 
mast. 

An  impressive  amount  of  space  has  been  devoted  to  scientific  purposes.  Five 
laboratories  and  offices  as  well  as  space  for  accomodating  portable  scientific  vans 
have  been  included  in  the  design.  There  is  berthing  for  10  scientists  and  techni¬ 
cians.  The  pride  of  the  oceanographic  instrumentation  is  the  computer  complex 
located  in  the  ocean  data  center,  which  assures  rapid  reduction  and  assimilation  of 
gathered  data  without  the  traditional  delay  until  return  to  homeport. 

The  aviation  capability  is  immediately  apparent  as  the  flight  deck  can  land 
the  heavy  HH3F  helicopter,  but  it  is  expected  that  two  HH52A  helicopters  will  make 
normal  deployments  with  POLAR  STAR.  The  helicopters  are  a  vital  tool  to  operating 
in  the  ice.  They  are  the  ship's  "eyes"  when  difficult  circumstances  are  encoun¬ 
tered. 

Perhaps  the  most  notable  feature  of  the  USCGC  POLAR  STAR  is  her  habitability. 
Careful  consideration  was  given  to  the  needs  of  the  crew  to  insure  they  received 
the  best  possible  in  care  and  comfort.  Duty  aboard  an  icebreaker  is  long  and  can  be 
arduous,  especially  when  your  ship  is  going  to  spend  eight  months  a  year  deployed 
at  sea.  The  manning  requirements,  for  example,  have  been  reduced  to  approximately 
150,  primarily  through  automation  and  the  use  of  low  maintenance  materials.  That 
number  compares  to  about  190  for  the  older  WIND  Class  icebreakers.  The  resulting 
extra  space  has  been  put  to  good  use.  Single,  double  and  four  man  staterooms  are 
provided  for  the  entire  crew.  Sizable  lounges,  a  library,  a  gymnasium  and  a  soda 
fountain  have  been  incorporated  into  the  ship.  Bright  colors  and  modern  decor  con¬ 
trast  noticeably  with  the  drabness  of  traditional  color  schemes.  These  and  other 
points  all  contribute  to  the  "livability"  of  POLAR  STAR. 

II.  COMMENTS  ON  VARIOUS  METHODS  OF  EVALUATING  AND/OR  PREDICTING  ICEBREAKER  PERFOR¬ 
MANCE  BY  CDR  G.  P.  VANCE 


A.  Introduction 


The  current  literature  contains  many  references  to  the  prediction 
of  icebreaker  performance.  Some  of  the  work  is  dated,  some  of  the  work  is  unique 
and  some  of  the  work  has  been  vastly  improved  upon.  This  note  will  attempt  to  clar¬ 
ify  the  utility  of  each  reference  in  the  prediction  technique,  however,  time  does 
not  permit  presentation  of  each  technique  and  sample  calculations.  Examples  can  be 
found  in  the  references  sited. 

The  prediction  techniques  can  be  broken  down  into  several  categories. 
Category  one  would  contain  the  techniques  that  are  analytical  in  nature  and  depend 
on  knowing  only  the  ship  and  ice  parameters,  i.e.,  no  model  tests  or  full  scale 
tests  data  are  necessary.  Included  in  this  category  are  predictions  for  continuous 
resistance  in  homogeneous  ice,  continuous  resistance  in  slush  ice,  ramming  capabi¬ 
lity  and  extraction  capability. 

The  second  category  involves  prediction  techniques  where  model  or  full 
scale  data  are  available  for  regression  analysis  to  establish  coefficients  for  the 
prediction  equation.  A  discussiorrof  the  various  techniques  is  contained  in  the 
following  sections.  A  summary  is  contained  in  Table  1. 

B.  Category  I 

1.  Discussion 


Kasteljan  (1968)  formulated  an  equation  that  expressed  ice  resistance 
encountered  by  an  icebreaker.  Although  his  equation  has  some  coefficients  in  it 
that  have  to  be  determined  experimentally,  the  values  published  are  taken  as  suit¬ 
able  for  all  hull  forms.  This  is  a  source  of  some  error  in  his  technique.  His 
equation  takes  the  form: 


are  determined  experimentally  and  the  other  terms  are  explained  in  his  1968  paper. 


The  equation  is  cumbersome  and  is  only  accurate  at  low  speeds  (Vance, 
1974).  It  appears  to  be  outdated  by  the  advance  of  technology  and  should  only  be 
used  as  a  check.  Details  of  its  use  are  contained  in  Kasteljan 's  paper  (1968). 

White  (1969)  studied  the  downward  force  generated  at  the  bow  of  an 
icebreaker  and  the  effects  that  the  bow  shape  had  on  this  force.  The  assumption 
is  that  the  icebreaker's  capability  is  proportional  to  bow  shape  alone,  which  is 
not  totally  true.  The  performance  is  dependent  on  the  entire  shape  of  the  vessel. 
Thus,  this  shortcoming  permeates  all  of  White's  work.  His  equations  should  be 
used  with  extreme  care,  particularly  the  simplification  to  ice  thickness  broken  for 
the  application  of  a  certain  downward  force. 

Examination  of  White's  equations  show  one  limitation  to  the  approach, 
i.e.,  it  does  not  take  into  consideration  the  full  ship.  It  should  be  used  with 
extreme  care  as  a  check. 

Milano's  (1973)  work  does  take  into  consideration  all  the  ship's 
parameters  as  well  as  the  ice  characteristics.  His  solution  is  a  complex  algo¬ 
rithm  that  must  be  run  on  a  computer.  It  has  been  my  experience  that  the  program 
does  relfect  what  is  happening  for  most  icebreaking  vessels.  However,  odd  changes 
in  shape  must  be  analyzed  with  care.  The  current  indications  are  that  Milano's 
program  yields  results  that  are  10  to  15  percent  lower  than  full  scale  results. 
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Table  1 


b 


Summary  of  I/B  prediction  techniques 
CONTINUOUS  HOMOGENOUS  ,CE 


AUTHOR 

Kastel jan 

White 

Milano 

Lewis 

Edwards 

Enkvist 

Vance 

CONTINUOUS  SLUSH  ICE 
llilano 
RAMMING 

White 
Lewi  s 
EXTRACTION 
White 


YEAR  PUBLISHED 

1968 

1969 

1973 

1971 

1972 
1972 

1974 

1975 

1969 

1971 

1969 


REMARKS 

Too  linear,  low 
Not  too  accureate 
Low,  needs  comp. 

Good  at  low  speeds 
Too  linear 
Good  at  low  speeds 
Good  within  5  percent 

Unproved 

Very  approximate 
Improved  White 

Relative  only 
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2.  Slush 

Milano's  (1975)  program  for  slush  ice  is  the  only  predictive  tool 
available  for  such  work.  It  was  first  presented  in  April  1975  and  has  not  been 
field  tested,  i.e.,  used  enough  by  design  engineers  to  give  confidence  in  the  re¬ 
sults.  There  is  not  much  that  can  be  said  for  its  application  to  date;  it  must  be 
tested  to  gain  more  insight  into  its  strength  and  weaknesses. 

3.  Ramming 

White  (1969)  developed  a  technique  that  yields  an  approximation  of 
the  thickness  of  homogeneous  ice  that  can  be  broken  in  a  ramming  mode.  Again  his 
development  is  limited  to  the  bow,  which  is  more  realistic  in  this  case.  Lewis 
(1971)  has  taken  White's  approach  and  made  some  improvements,  particularly  in  the 
application  of  the  downward  force  necessary  to  break  the  ice.  The  equations  should 
be  taken  as  approximates  and  are  sensitive  to  the  friction  coefficient  and  bow 
angles. 


4.  Extraction 

White's  (1969)  approach  to  extraction  difficulty  is  a  relative  one  and 
should  not  be  used  as  an  absolute  indication  of  the  ship's  ability  to  extract  itself 
with  a  given  thrust.  It  is  very  sensitive  to  bow  angle  and  one  should  be  aware  that 
he  is  referring  to  the  bow  angle  at  the  end  of  impact  rather  than  the  bow  angle  at 
the  designer's  waterline.  A  detailed  explanation  is  available  in  the  report  on  ram¬ 
ming  and  extraction  prediction  by  Vance  (1974). 

C.  Category  II 

Continuous  homogenous  ice: 

The  comments  stated  earlier  pertaining  to  Kasteljan's  equation  apply 
here  also.  The  coefficients  in  his  equation  are  found  through  model  tests  or  from 
full  scale  test  data  regression.  His  work  although  good  in  its  basic  concept  is 
not  precise  enough  for  today's  technology. 

Lewis,  in  1971,  proposed  an  equation  that  appeared  to  be  a  good  predic¬ 
tion  but  did  not  fit  all  shaped  ships.  It  depended  on  model  data  to  obtain  coeffi¬ 
cients  in  the  equation. 

Its  difficulty  lies  in  the  fact  that  it  does  not  contain  the  ship's 
length  and  the  beam  is  not  reflected  in  the  breaking  term.  Lewis,  et  al . ,  abandoned 
the  equation  in  1972  in  favor  of  another  equation  (Edwards,  1972). 

This  equation  has  the  inherent  weakness  of  being  linear  in  velocity 
which  is  contrary  to  most  field  data.  It  is  a  fairly  good  approximation  at  low 
speeds  where  the  velocity  effect  is  not  predominant. 

Enkvist  (1972)  used  another  regression  technique  that  he  felt  was  appli¬ 
cable  to  prediction  of  ice  resistance.  He  did  not  include  the  length  effect  and 
his  friction  effect  was  ambiguous.  His  technique  is  complicated  and  requires  several 
series  of  model  tests  to  determine  the  applicable  coefficients.  His  results  appear 
to  have  the  wrong  slope  for  the  resistance  curve  (Vance,  1974).  In  summary,  I  would 
not  recommend  Enkvist  technique  or  equations  for  predicting  full  scale  data. 

Vance  (1974)  proposed  a  regression  technique  that  included  the  length 
effect;  however,  because  of  the  lack  of  any  data  to  the  contrary,  he  assumed  the 
pressure  effect  was  neutral,  i.e.,  no  adverse  pressure  fields  existed.  He  also  as¬ 
sumed  that  the  friction  effect  was  spread  throughout  all  components  and  could  not  be 
separated  out.  His  results  were  within  5  percent  of  full  scale  data. 

It  must  be  kept  in  mind  that  none  of  these  techniques  include  the 
effect  of  pressure  and  are  somewhat  blase  about  the  effect  of  friction,  however. 
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that  is  the  state  of  the  art  to  date. 

D.  Summary 

In  summary,  I  would  recommend  Milano's  approach  for  vessels  that  have  not 
been  model  tested  or  are  still  on  the  drawing  board.  I  would  recommend  Vance's 
approach  for  analyzing  all  model  data. 


III.  PREDICTED  OPEN  WATER  PERFORMANCE  CHARACTERISTICS  OF  THE  POLAR  STAR 
BY  CDR  G.  P.  VANCE 

A.  Abstract 

The  enclosed  information  will  provide  an  estimate  of  the  performance  of 
the  POLAR  STAR  in  open  water.  The  information  is  based  on  the  Model  Test  Report 
P-223-H-09,  "Powering  Predictions  and  Flow  Observations  for  the  U.S.  Coast  Guard 
Proposed  Icebreaker  Design  (M-14-3)  Represented  by  Model  5245,"  by  M.P.  Lasky,  L.R. 
Crook  and  P.B.  Mathis,  dated  April  1971. 

The  RPM  reported  in  the  report  may  not  be  the  RPM  experienced  in  the  field, 
therefore  one  should  use  care  in  computing  the  torque,  i.e.,  each  torque  computation 
should  be  made  using  measured  RPM. 

One  should  also  be  careful  in  determing  the  power  distribution  when  under¬ 
way.  The  power  distributed  to  the  shafts  is  up  to  the  operator  and  may  not  be  equal, 
however,  the  total  power  should  be  equal  to  the  figures  presented  herein.  (See 
Table  2  and  Figs.  1-2). 
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B.  Basic  Equations 
PI  AN 

IHP - —  in  engine 

BHP  =  -  at  engine  coupling 


at  shaft  before  prop 

prop  (takes  into  account  stern  tube  losses) 

Thp  =  X  ^ 

550  powered  by  props 
RV 

EHP  =  g0g—  power  needed  to  drive  ship 

P.E.  -  Propulsive  Eff.  =  EH[P  (for  diesels) 

'IHP 

EHP 

=  -ggp  (for  turbines) 

Reference  page  372  PNA 

Propulsive  Coefficient  =  EHP  _  EHP  x  THP 

(or  Quasi-Propul sive  Coefficient)  PHP  ”  THP  PHP 


SHP  = 


d4G6M 

613033L, 


DHP  or  PHP  =  at 
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POLAR  STAR 
OPEN  WATER 
RESIST 


FIGURE  1  VELOCITY  (FEET  PER  SECOND) 


ouu 


; 


n 


VELOCITY  ( KTS) 
Figure  2. 


=  dh  x  ep  x  err 


Propulsive  Efficiency  -  EHP  PHP 

PHP  SHP 

=  eh  x  ep  x  err  x  et 

et  '  98%  mach.  aft  (rough  approximation) 

97%  mach.  midships  (rough  approximation) 

SHP  should  be  measured  as  close  to  stern  tubes  as  possible 

V-Va  =  wake  speed  v_Va 

=  wake  fraction  =  — y - 

Va  =  v ( 1-w ) 
l+io  is  wake  factor 

r  =  (l  -  t)T  total  resistance  or  available  thrust 
t  =  T"Rt  =  1  -  RT  thrust  deduction  fraction 

IT 

(1-t)  =  thrust  deduction  factor 


EHP  =  RtV  =  l-t_  =  e  =  hull  eff. 

THP  j -  1-w  h 

VA 


i.e.  want  t  small 


large 


VA  open  water  prop  eff. 
2ttNQ. 


eL  =  T  VA  eff.  of  prop  behind  ship 
b  2 NQ 


=  eK  =  Q, 


o  rel .  rotative  eff. 


.95<er  1.1  (rough  approximation) 

Dr  _  nor  -  EHP  _  RtV  _  RtV 
PC  or  QPC  dud  ~  I _ 


To  go  from  shaft  thrust  to  ship  resistance  we  need 


•  hull  efficient,  i.e.,  1-t  i.e.,  need  t  and  w 

1-0) 

•  need  ep,  i.e.,  open  water  prop  eff. 

•  need  er,  i.e.,  relative  rotative  eff. 

•  need  e^,  i.e.,  transmission  eff. 

C.  Expected  values  of  thrust 

The  total  thrust  that  can  be  expected  at  the  shafts  is  proportional  to 
the  total  resistance  encountered  by  the  ship.  If  we  assume  the  propulsion  coeffi¬ 
cient  presented  in  report  P-223-H-09  is  correct  then  the  expected  total  thrust  is 
presented  in  the  enclosed  table. 

It  should  be  emphasized  that  the  thrust  figure  is  the  total  and  the  thrust 
on  each  shaft  will  depend  on  the  power  distribution  set  by  the  operator,  i.e.,  1/3, 
1/3,  1/3  or  1/2,  0,1/2,  etc. 

Example  calculation 

Thrust  (at  shaft)  =  RESISTANCE 

P.C. 


RESISTANCE  =  EHPX550 
V 

R  (16  Kts )  =  5200X550  =  105,769  lbs. 

27.04 

Thrust  =  105769  =  176282  lbs. 

(Total)  .6 

D.  Expected  values  of  shaft  torque 
If  angle  of  twist  is  measured: 

shp  =  (<g4|5eN 

613033  Ls 

G-shear  modules  can  be  taken  as  11,900,000  PSI  for  steel  shafts 
N  is  revolutions  per  second 

Ls  length  over  which  angle  9  is  measured  in  inches 
d$  shaft  diameter  in  inches 
o  angle  of  twist  in  degrees  over  L$ 


For  torque  (foot  pounds) 


Q  in  ft  lbs 
N  in  RPS 


RPS  =  RPM 
60 


NOTE:  The  same  cautions  exist  for  torque  as  for  thrust,  i.e.,  the 
power  distribution  and  RPM  must  be  known. 

E.  Example  of  torque  calculation 

Q  =  SHP  550  _  SHP  87.535 
2ttN  N 


Assuming  RPM  of  175,  ..e.,  2.92  RPS  and  calculating  total  torque 


0  =  SHP  87.535 
2.92 


29.98  SHP 


At  8  KTS  open  •*  icr 

Q  =  1000  (29.90)  «  29977.74  «  2978  ft.  lbs. 

F.  Propeller  performance 

Currently  no  detailed  (i.e.,  computer  printout)  information  is  available 
for  the  controllable-pitch  propellers  on  the  POLAR  STAR.  The  enclosed  tables  and 
graph  will  give  some  indication  of  the  degradation  of  propeller  performance  with 
decreasing  speed  in  the  tow  rope  pull  condition.  One  must  be  careful  in  approxi¬ 
mating  the  P.C.  in  these  conditions,  particularly  in  ice.  Not  only  will  t  and  w 

change  with  speed  but  they  also  will  be  effected  by  ice  flow  in  a  detrimental 
manner.  The  extent  of  the  effect  is  difficult  to  estimate  due  to  the  lack  of  any 

model  or  full  scale  data.  (See  Fig.  3  and  Tables  3-7). 


IV.  PREDICTED  ICEBREAKING  PERFORMANCE  CHARACTERISTICS  OF  THE  POLAR  STAR 
BY  CDR  G.  P.  VANCE 


A.  Abstract 


The  curves  of  ice  resistance  versus  velocity  for  the  POLAR  STAR  have  been 
derived  utilizing  the  basic  program  of  Dr.  V.R.  Milano  (Resistance  to  Ship  Motion 
in  Sheet  Ice).  Caution  should  be  used  when  evaluating  changes  in  resistance  when 
there  has  been  a  dramatic  change  in  any  environmental  parameter.  The  user  should 
read  "Variation  of  Ship/Ice  Parameters  on  Ship  Resistance  to  Continuous  Motion  in 
Ice"  by  V.R.  Milano,  (April  1975),  SNAME  Eastern  Canadian  Section  paper,  before 
using  the  program.  The  program  has  been  modified  to  run  on  the  U.S.  Naval  Under¬ 
water  Systems  Center  UNIVAC  1108  computer  and  is  on  file  there.  The  program  can  be 
called  from  the  CG  R&DC  terminal  with  output  at  the  CG  R&DC  terminal  (very  slow) 
and/or  output  at  NUSC. 

The  program  plots  resistance  versus  velocity  curves  and  thrust  versus 
velocity  curves  (Fig.  4)  for  ice  thickness  from  1  to  9  feet  and  ship  speeds  from 
0  to  24  feet  per  second.  The  environmental  constants,  with  the  exception  of  flex¬ 
ural  ice  strength,  compressive  ice  strength,  the  coefficient  of  friction  and  the 
ice  specific  gravity  are  in  the  main  program  and  are  set  for  sea  ice.  The  program 
will  ask  for  the  values  of  flexural  ice  strength  (c,),  in  pounds  per  square  foot 
(PSF),  compressive  ice  strength  (o  )  in  PSF,  the  coefficient  of  friction  (u)  and 
specific  gravity  of  sea  ice.  Exmamples  of  typical  runs  are  contained  in  Table  8. 
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STOCK  PROPELLER 
PERFORMANCE 
OH 

POLAR  STAR 


FREE  ROUTE 


CENTER 


VELOCITY  (KNOTS) 


Table  3 

The  Test  Program  for  the  United  States  Coast  Guard  Proposed 
Icebreaker  Design  (M-14-3) 


Test 

Number 

Test 

Type 

Draft 

Figure 

Number 

Remarks 

1 

LOF 

28' 

5 

Bossing  Orientation 

•  . 

2 

EHP 

28' 

6  ,  7 

Bare  Hull 

m 

3 

EHP 

31' 

10,  19 

Free-Route 

m 

4 

SHP 

3 1 

10 

Free-Route,  Equal  SHP/Shaft 

■ 

5 

EHP 

28' 

11,  12,  20 

Free-Route 

m 

6 

SHP 

28' 

11,  12 

Free-Route,  Equal  SHP/Shaft 

7 

TRP-3 

28' 

13 

Equal  SHP  to  All  Shafts 

8 

TRP-3 

28' 

14 

No  SHP  to  Center  Shaft 

4 

. 

9 

TRP-3 

28' 

15 

No  SHP  to  Wing  Shafts 

10 

TRP-6 

28' 

13 

Equal  SHP  to  All  Shafts 

11 

TRP-6 

28' 

14 

No  SHP  to  Center  Shaft 

* 

12 

TRP-6 

28' 

15 

No  SHP  to  Wing  Shafts 

13 

TRP-9 

28' 

13 

Equal  SHP  to  All  Shafts 

14 

TRP-9 

28' 

14 

No  SHP  to  Center  Shaft 

S 

15 

TRP-9 

28' 

15 

No  SHP  to  Wing  Shafts 

16 

BP 

28' 

13 

Ahead,  Equal  SHP  to  All  Shafts 

17 

BP 

28' 

14 

Ahead,  No  SHP  to  Center  Shaft 

18 

BP 

28' 

15 

Ahead,  No  SHP  to  Wing  Shafts 

19 

BPA 

28' 

16 

Astern,  Equal  SHP  to  All  Shafts 

■ 

20 

BPA 

28’ 

17 

Astern,  No  SHP  to  Center  Shaft 

• 

21 

BPA 

28' 

18 

Astern,  No  SHP  to  Wing  Shafts 

22 

CWC 

28  7  31’ 

19-37 

Flow  Observations  of  Tests  3  through  21 

« 
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Table  4 


A  brief  description  of  the  stock  propellers  used  in  the  Testing 
Program  for  Simulating  Ahead  Modes  of  Operation  for  the  United 
statps  Coast  Guard  Proposed  Icebreaker  Design  (M-14-3; 


Propeller 

Number 

Shaft 

4305 

Starboard 

Wing 

4306 

Port 

Wing 

4307 

Center 

4454 

Starboard 

Wing 

4455 

Port 

Wing 

4456 

Center 

Shaft  Direction 
of  Rotation 

Pitch 

Diameter 

Outboard 

Right  Hand 

0.80 

Outboard 

Left  Hand 

0.80 

Right  Hand 

0.75 

Outboard 

Right  Hand 

1.00 

Outboard 

Left  Hand 

1.00 

Right  Hand 

1.00 

t 


I 


I 


Table  5 

Predicted  Free-Route  Propulsion  Coefficients 

(Test  4)  Draft  31  feet  (Zero  Trim)  Displacement  =  12,650  Tons 
Propellers  4305,  4306,  4307 


Vs 

KtS 

EHP 

SHP 

t 

C 

WT 

W 

c  WQ 

C 

JT 

W 

c 

W 

Thrust 

C 

(lbs) 

U 

16 

.60 

.  15 

.17 

.08 

.17 

.11 

.66 

.72 

.62 

.67 

48,300 

83,300 

17 

.60 

.15 

.17 

.08 

.17 

.10 

.65 

.71 

.62 

.67 

59,200 

102,000 

10 

.59 

.16 

.17 

.07 

.17 

.09 

.64 

.70 

.62 

.67 

76,800 

132,600 

19 

.58 

.17 

.17 

.05 

.17 

.08 

.62 

.69 

.62 

.67 

97,900 

169,000 

(Test  6) 

Draft  28  feet  (Zero  Trim) 
Propellers  4454,4455, 

Displacement  = 
4456 

10,870 

Ton 

Vs 

KtS 

EHP 

SHP 

t 

C 

WT 

W 

c  WQ 

W 

C 

JT 

W 

c  eP 

W 

Thrust 

C 

(lbs) 

W  ' 

10 

.65 

.14 

.17 

.06 

.18 

.08 

.83 

.89 

.73 

.70 

16,200 

27,600 

12 

.65 

.14 

.17 

.06 

.18 

.08 

.83 

.89 

.73 

.70 

24,000 

40,800 

14 

.65 

.  14 

.17 

.06 

.18 

.08 

.83 

.89 

.73 

.70 

33,800 

57,400 

15 

.65 

.14 

.17 

.06 

.18 

.08 

.83 

.89 

.73 

.70 

38,900 

66,400 

16 

.65 

.  15 

.17 

.06 

.19 

.08 

.83 

.89 

.73 

.70 

45,600 

78,100 

17 

.64 

.15 

.17 

.06 

.19 

.08 

.82 

.87 

.73 

.70 

56,300 

97,000 

IS 

.63 

.16 

.17 

.06 

.19 

.08 

.79 

.85 

.72 

.70 

73,400 

126,900 

14 

.61 

.17 

.17 

.05 

.19 

.08 

.76 

.83 

.71 

.69 

92,200 

159,100 

(Test  6A) 

Draft  28  feet  (Zero  Trim)  Displacement  = 
Propellers  4305,  4306,  4307 

10,870  Tons 

Vs 

l'ts 

EHP 

SHP 

t 

C 

WT 

w 

c 

W 

C 

JT 

W 

c 

U 

Thrust 

C 

(lbs) 

W 

K) 

.60 

.  Id 

.  17 

.08 

.16 

.10 

.67 

.74 

.62 

.67 

16,200 

28,000 

12 

.60 

.15 

.  17 

.08 

.16 

.10 

.67 

.74 

.62 

.67 

24,000 

41,400 

14 

.60 

.15 

.  17 

.08 

.16 

.10 

.67 

.74 

.62 

.67 

33,800 

58,200 

15 

.60 

.  15 

.17 

.08 

.17 

.10 

.67 

.74 

.62 

.67 

38,900 

67,000 

16 

.  6U 

.  15 

.17 

.08 

.17 

.10 

.67 

.73 

.62 

.67 

45,600 

78.500 

17 

.60 

.  15 

.  17 

.08 

.  18 

.10 

.66 

.72 

.62 

.67 

56,300 

97,000 

10 

.59 

.  16 

.  17 

.07 

.18 

.09 

.64 

.71 

.62 

.67 

73,400 

126,900 

19 

.53 

.  I7 

.17 

.05 

.18 

.07 

.62 

.70 

.62 

.67 

92,200 

159,100 

r 


i 


r 


r 


i 


t 
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Table  6 


Predicted  Tow-Rope  Pull  Propulsion  Coefficients 
Draft  =  28  feet  (Zero  Trim)  Displacement  =  10,870  Tons 
Equal  Power  Ahead  to  all  Shafts  -  Vs  =  3  Knots  -  Test 


TRP  (lbs) 

Thrust 

C 

(lbs) 

W 

R/T 

WT 

WQ 

JT 

C 

W 

C 

W 

C 

100,000 

41,100 

69,900 

.93 

.30 

-.10 

.15 

.02 

.21 

200,000 

78,200 

140,300 

.93 

.23 

-.18 

.09 

.00 

.17 

300,000 

108,300 

217,600 

.93 

.18 

-.26 

.05 

-.02 

.15 

400,000 

152,700 

280,800 

.93 

.14 

-.30 

.04 

-.02 

14 

500,000 

189,500 

351,600 

.93 

.12 

-.32 

.04 

.00 

.13 

600,000 

226,900 

421,700 

.93 

.10 

-.33 

.04 

.04 

.12 

700,000 

264,000 

492,100 

.93 

.06 

-.33 

.02 

.09 

.12 

800,000 

313,100 

550,600 

.93 

.03 

-.33 

.02 

.14 

.12 

876,000 

329,100 

616,400 

.93 

.01 

-.33 

.01 

.18 

.11 

Equal 

Power  Ahead 

to  All  Shafts 

-  Vs 

=  6  Knots 

-Test 

TRP  (lbs) 

Thrust 

C 

(lbs) 

W 

R/T 

WT 

WQ 

J. 

C 

W 

C 

w 

C 

100,000 

45,000 

79,800 

.91 

.24 

-.09 

.24 

.07 

.38 

200,000 

84,400 

150,100 

.91 

.26 

-.04 

.24 

.07 

.29 

300,000 

123,500 

220,600 

.91 

.29 

-.05 

.24 

.07 

.24 

400,000 

162,500 

291,300 

.91 

.30 

-.09 

.24 

.06 

.21 

500,000 

200,900 

362,500 

.91 

.30 

-.11 

.24 

.06 

.19 

600,000 

239,300 

433,800 

.91 

.29 

-.12 

.24 

.07 

.18 

700,000 

277,400 

505,300 

.91 

.27 

-.12 

.24 

.08 

.17 

800,000 

315,600 

576,800 

.91 

.26 

-.12 

.24 

.09 

.17 

810,000 

319,500 

583,900 

.91 

.26 

-.12 

.24 

.09 

.16 

17 
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Table  6  (Continued) 


f 


< 


< 


< 


1 


i 


Equal 

Power  Ahead  to 

all 

Shafts 

-  Vs  = 

9  1 

<nots  - 

Test 

13 

TRP  (lbs) 

Thrust 

C 

(lbs) 

W 

R/T 

W 

T 

WQ 

JT 

e 

P 

r 

C 

U 

C 

W 

C 

W 

C 

w 

100,000 

54,600 

89,800 

.90 

.24 

.03 

.19 

.09 

.46 

.58 

.66 

.60 

200,000 

93,500 

162,100 

.90 

.25 

.04 

.20 

.09 

.38 

.49 

.62 

.56 

r 

300,000 

132,300 

234,400 

.90 

.25 

.02 

.21 

.09 

.33 

.44 

.56 

.51 

400,000 

170,900 

306,800 

.90 

.26 

.00 

.24 

.09 

.30 

.40 

.51 

.47 

500,000 

209,500 

379,400 

.90 

.27 

.01 

.25 

.10 

.28 

.38 

.48 

.43 

r 

600,000 

248,100 

452,000 

.90 

.27 

.02 

.25 

.10 

.26 

.36 

.45 

.41 

700,000 

286.600 

524,500 

.90 

.28 

.03 

.25 

.10 

.24 

.34 

.43 

.39 

742.000 

304,700 

553,200 

.90 

.28 

.03 

.25 

.10 

.23 

.33 

.42 

.38 

Equal 

Power  Ahead  to 

All 

Shafts 

-  Vs  = 

0  Knots  - 

Test 

16 

TRP  (lbs) 

Thrust 

C 

(lbs) 

W 

R/T 

WT 

WQ 

JT 

ep 

C 

W 

C 

W 

C 

W 

C 

100,000 

36,800 

67,100 

.96 

... 

... 

... 

•A  — 

200,000 

72,900 

133,700 

.96 

--- 

... 

... 

... 

300,000 

108,300 

202,200 

.96 

— 

... 

... 

... 

400,000 

142,400 

274,300 

.96 

— 

... 

... 

... 

500,000 

179,700 

340,300 

.96 

--- 

... 

... 

... 

600,000 

217,900 

405,100 

.96 

— 

— 

... 

... 

700,000 

255,600 

475,300 

.96 

— 

— 

... 

... 

000,000 

292,600 

547,400 

.96 

--- 

— 

... 

... 

900,000 

329,500 

618,000 

.96 

— 

... 

... 

... 

926,000 

338,600 

635,100 

.96 

— 

... 

— 

— 

i 
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Table  6  (Continued) 


Equal  Power 

Ahead  to  Wing  Shafts  Only 

-  Vs  = 

3  Knots 

-  Test  8 

TRP  (lbs) 

Thrust  (lbs) 

R/T 

WT 

WQ 

JT 

6P 

* 

100,000 

115,200 

.93 

-.09 

.05 

.27 

.35 

200,000 

223,400 

.93 

-.05 

.12 

.20 

.26 

300.000 

331,300 

.93 

-.05 

.18 

.16 

.21 

r 

400,000 

439,500 

.93 

-.07 

.20 

.15 

.18 

500,000 

543,700 

.93 

-.07 

.26 

.13 

.17 

577,500 

631,400 

.93 

-.08 

.32 

.13 

.16 

r 

Equal  SHP  to  Wing  Shafts  Only 

-  Vs  =  6 

Knots 

-  Test 

11 

TRP  (lbs) 

Thrust  (lbs) 

R/T 

WT 

WQ 

JT 

6P 

■ 

r 

100,000 

130,000 

.93 

-.01 

.07 

.43 

.52 

• 

200,000 

237,400 

.93 

-.03 

.10 

.34 

.42 

300,000 

344,100 

.93 

-.04 

.11 

.30 

.37 

» 

400,000 

451,500 

.93 

-.06 

.12 

.27 

.34 

500,000 

558,200 

.93 

-.08 

.14 

.25 

.32 

529,000 

589,300 

.93 

-.09 

.16 

.24 

.31 

l 

Equal  Power 

Ahead  to  Wing  Shafts  Only 

-  Vs  = 

9  Knots 

Test  14 

TRP  (lbs) 

Thrust  (lbs) 

R/T 

WT 

WQ 

JT 

6P 

1 

100,000 

151,100 

.93 

.03 

.10 

.50 

.58 

1 

200,000 

259,700 

.93 

.02 

.10 

.42 

.51 

300,000 

367,500 

.93 

.00 

.11 

.38 

.46 

w  ' 

400,000 

476,000 

.93 

-.02 

.12 

.35 

.43 

468,800 

495,800 

.93 

-.02 

.12 

.34 

.43 

t 
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Table  6  (Continued) 

Equal  Power  Ahead 

to  Wing  Shafts  Only 

-  Vs 

=  0  Knots 

-  Test  17 

TRP  (lbs) 

Thrust  91 bs ) 

R/T 

WT 

“q 

JT 

100,000 

103,100 

.97 

— 

— 

— 

200,000 

206,200 

.97 

— 

— 

— 

300,000 

309,300 

.97 

— 

— 

— 

400,000 

412,400 

.97 

— 

— 

— 

500,000 

515,500 

.97 

— 

— 

— 

600,000 

618,600 

.97 

— 

— 

— 

620,000 

639,200 

.97 

— 

— 

— 

Ahead  Power  to 

Center  Shaft  Only  - 

Vs  = 

3  Knots  - 

Test  9 

TRP  (lbs) 

Thrust  (lbs) 

R/T 

WT 

"q 

JT 

50,000 

66,600 

.92 

.25 

.22 

.18 

100,000 

121,300 

.92 

.25 

.22 

.14 

150,000 

175,900 

.92 

.24 

.21 

.12 

200,000 

229,300 

.92 

.24 

.20 

.11 

250,000 

281,800 

.93 

.23 

.18 

.10 

300,000 

334,200 

.93 

.22 

.17 

.10 

301,000 

335,400 

.93 

.22 

.17 

.10 

Ahead  Power  to 

Center  Shaft  Only  - 

Vs  = 

6  Knots  - 

Test  12 

TRP  (lbs) 

Thrust  (lbs) 

R/T 

WT 

WQ 

JT 

50,000 

87,200 

.91 

.26 

.16 

.28 

100,000 

142,400 

.91 

.27 

.17 

.23 

150,000 

197,300 

.91 

.27 

.18 

.20 

200,000 

250,700 

.91 

.27 

.19 

.18 

250,000 

303,800 

.92 

.27 

.19 

.17 

267,000 

322,300 

.92 

.27 

.19 

.16 

E 


Table  6  (Continued) 


Ahead  Power 

to  Center  Shaft  Only  - 

Vs  =  9 

Knots  - 

Test  15 

TRP  (lbs) 

Thrust  (lbs) 

R/T 

WT 

WQ 

JT 

6P 

50,000 

128,600 

.79 

.14 

.08 

.39 

.48 

100,000 

181,900 

.83 

.13 

.07 

.34 

.42 

150,000 

236,700 

.85 

.13 

.06 

.31 

.38 

200,000 

290,800 

.86 

.13 

.04 

.29 

.36 

213,000 

305,200 

.86 

.13 

.04 

.28 

.35 

Ahead  Power 

to  Center  Shaft  Only  - 

Vs  =  0 

Knots  - 

Test  18 

TRP  (lbs) 

Thrust  (lbs) 

R/T 

WT 

WQ 

JT 

6P 

50,000 

51,500 

.97 

— 

— 

— 

— 

100,000 

103,100 

.97 

— 

— 

— 

— 

150,000 

154,600 

.97 

— 

— 

— 

— 

200,000 

206,200 

.97 

— 

— 

— 

— 

250,000 

257,700 

.97 

— 

— 

— 

— 

300,000 

309,300 

.97 

— 

— 

— 

— 

337,000 

347,400 

.97 

— 

— 

— 

— 

Equal  Power  Astern  to 

All  Shafts 

-  Vs  =0 

Knots  - 

Test  19 

TRP  (lbs) 

Thrust 

C 

(lbs) 

W 

R/T 

WT 

WQ 

JT 

6P 

100,000 

38,100 

70,100 

.93 

— 

— 

— 

— 

200,000 

78,800 

146,100 

.89 

— 

— 

— 

— 

300,000 

119,900 

226,400 

.87 

— 

— 

— 

— 

400,000 

159,600 

306,300 

.86 

— 

— 

— 

— 

500,000 

200,000 

381,800 

.86 

— 

— 

— 

— 

600,000 

238,000 

459,700 

.86 

— 

— 

— 

— 

700,000 

276,300 

536,400 

.86 

— 

— 

— 

— 

800,000 

315,300 

613,600 

.86 

— 

— 

— 

— 

892,000 

352,700 

683,500 

.87 

— 

— 

— 

— 

21 
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1  Table  6  (Continued)  j 

Equal  Power  Astern  to  Wing  Shafts  Only  -  Vs  =  0  Knots  -  Test  20 


TRP  (lbs) 

Thrust  (lbs) 

R/T 

WT 

WQ 

JT 

6P 

\ 

100,000 

119,000 

.84 

— 

— 

— 

— 

200,000 

240,600 

.83 

— 

— 

— 

— 

300,000 

362,500 

.83 

— 

— 

— 

— 

! 

400,000 

484,400 

.83 

— 

— 

— 

— 

: 

500,000 

606,000 

.83 

— 

— 

— 

— 

1 

561,000 

686,000  .83 

Astern  Power  to  Center  Shaft  Only  - 

Vs  =  0 

Knots  - 

Test  21 

~  ” 

TRP  (lbs) 

Thrust  (lbs) 

R/T 

WT 

WQ 

JT 

6P 

50,000 

61,600 

.81 

— 

— 

— 

— 

1 

100,000 

124,800 

.80 

— 

— 

— 

— 

150,000 

187,100 

.80 

— 

— 

— 

— 

200,000 

247,300 

.81 

— 

— 

— 

— 

\ 

250,000 

307,400 

.81 

— 

— 

— 

— 

290,000 

355,700 

.81 

— 

— 

— 

— 

i 
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ft 


Table  7 


Summary  of  Predicted  Towing  Performance  for  the  M-14-3  Design  Icebreaker 
Shaft  Horsepower  =  20,000  SHP/SHAFT 


Maximum 

Power  to 

All  Shafts 

VS 

RPM 

TRP  (lbs) 

Thrust 

(lbs) 

Total  SHP 

KTS 

C 

W 

C 

W 

0  Astern 

164.5 

157.8 

892,000 

352,700 

683,500 

60,000 

0  Ahead 

165.6 

158.6 

926,000 

338,600 

635,100 

60,000 

r 

3  Ahead 

167.1 

163.7 

876,000 

329,100 

616,400 

60,000 

6  Ahead 

170.0 

168.3 

810,000 

319,500 

583,900 

60,000 

9  Ahead 

176.9 

175.8 

742,000 

304,700 

553,200 

60,000 

r 

Maximum  1 

Power  to  Wing  Shafts  - 

No  Power  to 

Center  Shaft 

0  Astern 

Wind  Milling 

158.9 

561,000 

686,200 

40,000 

0  Ahead 

Wind  Milling 

157.8 

620,000 

639,200 

40,000 

r 

3  Ahead 

Wind  Milling 

162.3 

577.500 

631,400 

40,000 

6  Ahead 

Wind  Milling 

168.8 

529,000 

589,300 

40,000 

9  Ahead 

Wind  Mil  lino 

175.2 

468,800 

547  ,400 

40 ,000 

i 

Maximum 

Power  to  Center  Shaft  - 

No  Power  to  Wing  Shafts 

0  Astern 

164.2  Wind 

Mill ing 

290,000 

355,700 

20,000 

0  Ahead 

162.7  Wind 

Milling 

337,000 

347,400 

20,000 

9 

3  Ahead 

166.8  Wind 

Milling 

301,000 

335,400 

20,000 

6  Ahead 

171.1  Wind 

Milling 

267,000 

322,300 

20,000 

9  Ahead 

176.9  Wind 

Milling 

213,000 

305,200 

20,000 

t 
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AR  TOW  ROPE  THRUST 


0.00  2.50  5.00  7.50  10.00  12.50  15. OC 

VELOCITY  IN  FPS 
THRUST  VS  VELOCITY 


Table  8 


Summary  of  runs 


Run  Ship 

af  (PSF) 

oc  (PSF) 

SG  ICE 

Remarks 

1.  POLAR  STAR 

10,000 

100,000 

.2 

1.03 

CG  R&DC  Ouput 

2.  POLAR  STAR 

10,000 

100,000 

.2 

3.  POLAR  STAR 

12,000 

100,000 

.2 

4.  POLAR  STAR 

15,000 

100,000 

.2 

5.  POLAR  STAR 

10,000 

100,000 

.4 

6.  POLAR  STAR 

10,000 

200,000 

.2 

7.  POLAR  STAR 

10,000 

100,000 

.2 

8.  POLAR  STAR 

10,000 

100,000 

.2 

Plot  RT  Change 
Plot  to  Milano 

9.  POLAR  STAR 

15,000 

100,000 

.2 

10.  POLAR  STAR 

20,000 

100,000 

.2 

11.  POLAR  STAR 

10,000 

100,000 

.2 

Plot  missing 

12.  WIND 

24,000 

40,000 

.2 

.92 

Same  as  Milano 

SG  ice  changed 
of  high,  ODR  RES 

13.  WIND 

10,000 

100,000 

.2 

14.  POLAR  STAR 

10,000 

100,000 

.2 

15.  POLAR  STAR 

20,000 

100,000 

.2 

16.  POLAR  STAR 

12,000 

100,000 

.2 

17.  POLAR  STAR 

15,000 

100,000 

.2 

18.  POLAR  STAR 

10,000 

100,000 

.1 

19.  POLAR  STAR 

15,000 

100,000 

.3 

Term  doubt 

20.  POLAR  STAR 

15,000 

100,000 

.3 

Term  doubt 

21.  POLAR  STAR 

15,000 

50,000 

.2 

Printout  missing 

22.  POLAR  STAR 

15,000 

50,000 

.3 

23.  WIND 

15,000 

50,000 

.2 

24.  MACKINAW 

15,000 

50,000 

.2 

25 


Table  8  (Continued) 


Run 

Ship 

of  (PSF) 

oc  (PSF) 

u 

SG  ICE 

Remarks 

25. 

MACKINAW 

20,000 

50,000 

.2 

.92 

\ 

26. 

MACKINAW 

30,000 

50,000 

.2 

.92 

27. 

POLAR  STAR 

30,000 

50,000 

.2 

.92 

28. 

POLAR  STAR 

15,000 

50,000 

.2 

.82 

CHD  pro  to  call 

1 

for  SPGR  ice 

29. 

POLAR  STAR 

15,000 

50,000 

.1 

.82 

30. 

POLAR  STAR 

20,000 

50,000 

.2 

.92 

1 

31. 

MACKINAW 

20,000 

50,000 

.2 

.92 

32. 

POLAR  STAR 

12,000 

50,000 

.2 

.92 

33. 

POLAR  STAR 

15,000 

50,000 

.2 

.96 

i 

34. 

POLAR  STAR 

15,000 

52,500 

.2 

.92 

35. 

POLAR  STAR 

15,000 

45,000 

.2 

.92 

36. 

POLAR  STAR 

14,300 

50,000 

.2 

.92 

\ 

37. 

POLAR  STAR 

16,000 

50,000 

.2 

.92 

38. 

POLAR  STAR 

18,000 

50,000 

.2 

.92 

39. 

POLAR  STAR 

15,000 

50,000 

.1 

.92 

m 

40. 

POLAR  STAR 

15,000 

50,000 

.25 

.92 

41. 

POLAR  STAR 

15,000 

50,000 

.235 

.92 

42. 

MACKINAW 

20,000 

50,000 

.2 

.945 

FR  water 

i 

43. 

POLAR  STAR  30’ 

15,000 

50,000 

.2 

.945 

FR  water 

44. 

POLAR  STAR  30' 

15,000 

50,000 

.2 

.92 

Salt  water 

45. 

POLAR  STAR  28' 

15,000 

50,000 

.2 

.92 

GML=3501 . 1 

G  BETA  14  ST 

i 

46. 

POLAR  STAR  30' 

15,000 

50,000 

.2 

.92 

G  DELT  14  ST 

« 
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B.  Introduction 


The  ship  data  is  read  into  the  program  as  data  file.  At  present  the  pro¬ 
gram  is  set  up  such  that  the  data  file  for  the  ship  of  interest  must  be  called  for. 
The  POLAR  STAR  data  is  contained  in  a  data  package  called  ICEDTA.  Any  ship's  data 
can  be  called  for  by  using  the  @  ASGA  (file  name)  and  the  @  USE  15,  (file  name), 
cards.  The  following  tables  describe  the  program  input  data,  with  examples  for 
the  POLAR  STAR  (Table  9),  the  WINDCLASS  (Table  10)  and  the  MACKINAW  (Table  11), 
the  program  constants  and  the  program  output.  The  wind  data  is  contained  in  a  file 
called  WINDTA,  the  MACKINAW  data  is  in  a  file  called  MACDTA. 

Once  the  terminal  is  activated  by  turning  it  on,  dialing  NUSC  (447-3251) 
and  entering  the  CG  R&DC  ID  U1108P,  the  runs  are  made  by  following  the  enclosed 
samples. 

The  printout  of  the  curves  of  Resistance  versus  Velocity  and  Thrust 
versus  Velocity  are  contained  in  appendix  A  of  this  report. 

C.  Nomenclature  for  ice  resistance  program  (Input) 


1.  RLEN 

Length  between  perpendiculars  in  feet 

2.  BX 

Maximum  beam  in  feet 

3.  DEL 

Displacement  in  tons 

4.  SHP 

Shaft  horsepower  per  shaft  in  horsepower 

5.  PDIA 

Propeller  diameter  in  feet 

6.  X 

Distance  from  bow  to  the  section  of  maximum  beam  in 
feet 

7 .  DELD 

Distance  from  maximum  beam  to  LCG  IN  FEET  (aft  is 
positive) 

8.  CBXL2 

Distance  from  maximum  beam  to  amidship  in  feet  (aft  is 
positive) 

9.  DRAFT 

Draft  in  feet 

10.  CX 

Midship  section  coefficient 

11.  CW 

Waterplane  coefficient 

12.  GML 

Longitudinal  metacentric  height  in  feet  at  stated  draft 

13.  ALPHA 

Angle  of  inclination  of  the  bow  measured  from  horizontal 
in  radians 

14.  ALCG 

Distance  from  LCG  to  LCF  in  feet  (aft  is  positive) 

15.  NP 

Number  of  propellers 

16-18  CBI  (I) 

Station  beam  coefficient,  C..=B./B  a  21  component 
array  from  station  0  to  20  01  1  x 

19-21.  CXI  (I) 

Station  area  coefficient  C  .=A  ./A  a  21  component  array 
from  station  0  to  20  xi  xi  x 

Table  9 


POLAR  STAR  input  data  (ICEDTA)  28  foot  draft 
♦TERMINAL  INACTIVE* 

>@RUN  TDSMET,  X832J0ESMITH ,A99614STC000,60,500 
DATE:  061775  TIME:  152842 
>@ASG,AX  ICEDTA. 

READY 

>@DATA,L  ICEDTA. 

1234567890123456789012345678901234567890123456789012345678901234567890123456789012 
25 :>  P  1  40 
352.0 
78.0 
10863.0 
20000.0 
16.0 
140.  g 
33.3 
35.2 
28.0 
.852 
.740 
350.1 
.26179 
1.35 

3 


0.0 

.242 

.465 

.648 

.789 

.890 

.955 

.989 

1.0 

.999 

.990 

.974 

.953 

.925 

.883 

.823 

.740 

.622 

.459 

.249 

.000 

0.0 

.032 

.150 

.353 

.556 

.735 

.856 

.922 

.973 

.992 

1.0 

.968 

.911 

.824 

.714 

.636 

.406 

.267 

.150 

.032 

.000 

.524 

.593 

.716 

.794 

.922 

1.114 

1.257 

1.295 

1.309 

1.326 

1.335 

1.344 

1.335 

1.274 

.489 

.471 

.419 

.353 

.257 

.175 

.112 

.056 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

1. 

2. 

3. 

4. 

5. 

6. 

7. 

8. 

9. 

0.0 

SCAN :4 
EOF : 35 
0:>3CALE  1 

1234567890123456789012345678901234567890123456789012345678901234567890123456789012 
0: >EXIT 

LINES : 35  FIELDATA 
>3FREE  ICEDTA. 

READY 

• @@SKIP  60 
-•OF  IN 
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Table  10 


WINDCLASS  input  data  (WINDTA) 

12345678901234567890123456789012345678901233456789012345678901234567890123456789012 

25  :>  P  1  40 

250.0 

63.75 
5300.0 
5000. 

17. 

125. 

1. 

0. 

25.75 
.752 
.724 
201.4 
.523 
-.220 

2 


0.0 

0.225 

.432 

.605 

.737 

.840 

.910 

958 

.982 

.992 

1.0 

.993 

.9770 

.951 

91 

.844 

.754 

.630 

.462 

.246 

0.0 

0.0 

.045 

.202 

.386 

.554 

.696 

.810 

.895 

.947 

.980 

1.0 

.997 

.978 

.929 

.854 

.741 

.593 

.415 

.250 

.079 

0. 

.916 

.916 

.960 

1.082 

1.169 

1.239 

1.265 

274 

1.274 

1.274 

1.274 

1.274 

1.274 

1.274 

.526 

.471 

.436 

.367 

.285 

.209 

.  148 

.087 

.052 

.026 
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0.0 

0.0 

0.0 
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Table  11 

MACKINAW  input  data  (MACDTA. ) 

1234567890123456789012345678901234567890123456789012345678901234567890123456789012 

25  :>  P  1  40 

280. 

70. 

5252. 

5000. 

17. 

140. 

1.0 

00. 

19. 

.812 
.728 
321.5 
.523 
.  1 

2 


0.000 

0.236 

0.434 

0.599 

0.732 

0.830 

0.906 

0.957 

0.983 

0.997 

1.000 

0.998 

0.984 

0.957 

0.917 

0.855 

0.769 

0.644 

0.473 

0.221 

0.0 

0.000 

0.052 

0.239 

0.411 

0.567 

0.699 

0.809 

0.891 

0.950 

0.985 

1.000 

0.993 

0.958 

0.906 

0.809 

0.676 

0.527 

0.351 

0.181 

0.047 

0.0 

0.843 

0.855 

0.890 

0.942 

0.994 

1.047 

1.094 

1.143 

1.178 

1.204 

1.213 

1.213 

1.204 

1.18 

0.532 

0.489 

0.436 

0.367 

0.300 

0.227 

0.166 

0.105 

0.061 

0.035 

0.018 

0.0 

0.0 

0.0 

1. 

2. 

3. 

4. 

5. 

6. 

7. 

8. 

9. 

0.0 

SCAN:  4 
EOF: 35 
0:  >SCALE  1 

1234567890123456789012345678901234567890123456789012345678901234567890123456789012 
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22-23.  GBETA  (I) 

24-25.  GDELT  (I 


GBETA  (I) 


26-34.  HICE 
35. 

D.  Nomenclature  for 

1.  RHOW 

2.  RHOI 

3.  GAMW 

4.  SPGR 

5.  SJG 

6.  SIGC 


7.  FRICT 

8.  XNU 

9.  GRA 

10.  DELCF 


Station  transverse  spread  angle  complement  in  radians 
an  11  or  15  component  array  station  0  to  10  or  14  taken 
from  body  plan  depending  on  location  of  maximum  beam 

Station  waterline  inclination  in  radian  an  11  or  14 
component  array  station  0  to  10  or  13  taken  from  half 
breadth  plan  depending  on  location  of  maximum  beam 

GDELT  (I) 


Ice  thickness  in  feet 

Termination  of  loop  with  0.0  ice  thickness 

ice  resistance  program  (Program  constants) 

Mass  density  of  water  in  LBsec2/FT4  normally  taken  as 
1.99 

2  4 

Mass  density  of  ice  in  LBsec  /FT  normally  taken  as  1.79 

3 

Specific  weight  of  water  normally  taken  as  65.0  LB/FT 
Specific  gravity  of  ice  normally  taken  .92 

O 

Tensile  (flexural)  strength  of  ice  in  LB/FT  .  In  the 
area  of  10,000  PSF  for  seawater  ice  and  20,000  PSF  for 
fresh  water  ice.  This  value  should  be  investigated  by 
user. 

2 

Compressive  strength  of  ice  in  LB/FT  .  This  value  de¬ 
pends  on  ice  temperature.  50,000  can  be  taken  as  an 
average  value. 

Coefficient  of  dynamic  friction.  This  value  will  vary 
and  should  be  investigated  by  user.  Range  of  values 
would  be  .1  to  .5. 

Kinematic  viscosity  of  water  normally  taken  as  1.97  x 
10"5  FTvsec 

2 

Acceleration  of  gravity  normally  taken  as  32.2  FT/sec 

Hull  roughness  allowance  normally  taken  as  .0004  (used 
for  open  water  resistance) 


11.  PSI0 


Bow  wedge  included  angle  normally  taken  as  1.18  radians 


i 


I 


1 


1 


12. 

PSIC 

Cusp  wedge  included  angle  normally  taken  as 

1.85  radians 

13. 

E 

Young  modulus  for  ice  normally  taken  as  1  x 

108  LB/FT2 

14. 

XK 

Foundation  (water)  modulus  normally  taken  as 

64  LB/FT2 

15. 

ACCYF 

An  epsilon  factor  used  in  computing  open  water  resis¬ 
tance  normally  taken  as  .001E-3 

Nomenclature  for 

ice  resistance  program  (Output) 

1. 

HICE 

Ice  thickness  in  feet 

2. 

U 

Ship  speed  in  feet  per  second 

3. 

El 

Energy  through  ice  filled  channel  in  foot  pounds 

4. 

THURST 

Total  thrust  available  in  pounds 

5. 

TT1 

Time  to  break 

6. 

E3 

Energy  for  climbing  on  ice  in  foot  pounds 

7. 

TT 

Time  in  seconds 

8. 

E4 

Energy  for  fracturing  ice  in  foot  pounds 

9. 

E5 

Energy  for  submerging  ice  in  foot  pounds 

10. 

E21 

Energy  for  local  crushing  in  foot  pounds 

11. 

ET 

Total  energy  in  foot  pounds 

12. 

FT 

Total  resistance  in  pounds 

13. 

XFT 

Non-dimensional  resistance  (R/pgBH2^ 

14. 

XFROU 

Thickness  froude  number  (U//gFi) 

15. 

XSIG 

Non-dimensional  strength  (o/pgh) 

POLAR  STAR  station  beam  coefficient 

STA 

BEAM/ 2 

Bi/Bx 

0 

0 

.000 

1 

9-5-2 

.242 

2 

18-1-5 

.465 

3 

25-3-1 

.648 

4 

30-9-0 

.789 

5 

34-8-2 

.890 
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;  • 

STA 

BEAM/ 2 

Bi/Bx 

6 

37-2-5 

.955 

7 

38-6-6 

.989 

T 

*8  (bx) 

38-11-7 

1.000  (38.9895) 

9 

38-11-6 

.999 

10 

38-7-2 

.990 

S 

11 

37-11-5 

.974 

12 

37-2-0 

.953 

13 

36-0-7 

.925 

1 

14 

34-5-0 

.883 

15 

32-1-0 

.823 

16 

28-10-3 

.740 

i 

17 

24-3-3 

.622 

18 

17-10-5 

.459 

19 

9-8-5 

.249 

20 

0  0  0 

.000 

G.  POLAR  STAR  station 

area  coefficient 

STA 

AREA 

Ai/Ax 

0 

0 

.000 

1 

60 

.032 

2 

280 

.150 

* 

3 

660 

.353 

4 

1040 

.556 

5 

1374 

.735 

i 

6 

1600 

.856 

7 

1725 

.922 

8 

1820 

.973 

■ ' 

9 

1855 

.992 
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SECTION  AREA  CURVE 


Figure 


STA 

14 

15 


ANGLE  (DEC) 
71.0 
66.8 


I.  Station  waterline  inclination  [GDELT  (I)] 


STA 


ANGLE  (DEG) 


0 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 
11 
12 

13 

14 

15 


28 

27 

24 

20.25 

14.70 

10.00 

6.4 

3.2 

0.0 

-1.0 

-2.0 

-2.5 

-3.5 

-5.0 

-7.0 

-10.0 


ANGLE  (RAD) 
1.239 
1.166 

ANGLE  (RADIANS) 

(Used)  (Act) 

.489 

.471 

.419 

.353 

.257 

.175 

.112 

.056 

.000 

.000  -.017 

.000  -.035 

.000  -.044 

-.061 
-.087 
-.122 
-.175 


U.  Expected  thrust  values  when  breaking  ice 

The  thrust  experienced  at  the  shafts  when  breaking  ice  will  be  directly 
proportioned  to  the  resistance  experienced  by  the  vessel.  The  expected  value  of 
total  thrust  (as  in  open  water,  the  distribution  of  power  on  the  shafts  is  up  to 
the  operator)  can  be  found  by  referring  to  the  resistance  curve  that  matches  the 
operating  conditions.  All  that  need  be  done  to  find  total  thrust  at  the  shaft  is 
divide  the  total  resistance  by  the  propulsive  coefficient. 

Example  calculation  from  Run  #21 
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Of  =  15,000  PSF  0C  =  50,000  PSF  u  =  .2  SG  =  .92 
R  =  340,000  lbs. 


c 


9 


I 


1 


i 


i 

t 


Tshaft  =  R  =  340,000 

K~  .6 

566,667  lbs.  (total ) 

NOTE:  The  P.C.  may  be  less  than  .6,  which  was  taken  from  the  open  water 
tests,  because  of  the  ice  interrupting  the  water  flow  about  the  propellers. 

K.  Expected  torque  values  when  breaking  ice 

When  in  the  icebreaking  mode,  the  torques  experienced  on  the  shafts  will 
be  much  higher  than  those  experienced  while  in  open  water.  The  torque  will  depend 
on  the  total  resistance  of  the  vessel  as  well  as  the  speed  through  the  ice. 

In  order  to  estimate  the  torque  the  power  distribution  (i.e.,  1/3,  1/3,  r 

1/3)  must  be  specified.  The  total  resistance  is  distributed  accordingly  (another 
assumption)  and  the  resistance  changed  to  EHP  which  is  converted  to  SHP,  using  a 
P.C.  of  approximately  .6  (another  assumption,  the  P.C.  will  probably  be  lower  be¬ 
cause  of  ice  interferring  with  water  flow).  The  SHP  can  be  converted  to  shaft 
torque  knowing  the  RPM  of  the  shaft. 

This  procedure  does  not  take  into  account  any  torques  that  may  be  en- 
countered  due  to  interaction  between  the  ice  and  the  propeller  itself.  These 
torques  (i.e.,  ice-propeller)  will  show  up  as  large  peaks  on  an  oscillograph  or 
recorder  record  and  can  be  eliminated  from  any  power  calculation  if  desired. 

L.  Calculation  of  expected  ice  resistance  torque 

r 

Q  =  550  SHP 
2  it  N 

N  -  RPS  S  -  FT  LBS 


P-C-  =  eh  ep  err  et 


EHP 


Ice  resistance  x  vel 
500 


Ice  resist  =  lbs  (obtained  from  Resist  Curves) 

VEL  -  FPS 

NOTE:  Ice  resistance  must  be  divided  in  accordance  with  power  distri¬ 
bution  example  calculation 

Ice  resistance  =  340,000  lbs 

From  plots  R  vs  V  for  run  #21 

13  FPS  (7.7  KTS)  4  ft  ice  of  =  15,000  PSF  u  =  .2 
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EH  P  = 

EPH 

3 


50,000  SG  =  .92  RPM  =  175 

340,000  x  13  _  fln,, 

550  8036 


8036 


^  2678.8 


SHP 


xs 


EHP 


2678.8 


=  4464.6 


P.C.  .6 
550  SHP  =  550  (4464.6) 


2-  N 


2tt  2.92 


Q  =  133906  ft  lbs  per  shaft  rotating  at  175  RPM 
M.  Prediction  program  verification 

Although  the  program  may  be  over  sensitive  to  the  various  parameters 
entered,  both  ship  and  environment,  one  can  obtain  an  appreciation  of  its  relative 
merit  by  comparing  the  output  for  the  POLAR  STAR  with  the  MACKINAW.  The  MACKINAW 
output  can  then  be  compared  to  full-scale  data  obtained  in  field  tests.  This  is 
done  for  Run  *21  and  Run  *  24  (i.e.,  thickness  plot  -  see  Figs.  6-7). 

The  comparison  must  be  made  with  care  since  Run  #21,  and  Run  *24  made 
with  the  following  parameters: 


7f  =  15,000  PSF 
:c  =  50,000  PSF 

■m  =  .  2 

SG  =  .92 


The  program  appears  to  predict  a  resistance  that  is  15  percent  lower  than 
full  scale  at  speeds  greater  than  7  FPS  (4.2  KTS).  There  may  be  several  reasons 
for  this  discrepancy: 

•  The  coefficient  of  friction  may  not  have  been  measured  correctly  in 

the  field 


•  The  environmental  or  power  parameters  may  not  have  been  correctly 
determined  in  the  field 

•  The  program  algorithms  may  not  be  correct. 

In  any  event,  one  should  be  aware  of  these  shortcomings  when  comparing 
field  data  with  the  program  output. 

V.  PREDICTION  OF  ERROR  SOURCES  FOR  FULL  SCALE  TEST  OF  POLAR  STAR  BY  G.  P.  VANCE 
A.  Abstract 

An  error  analysis,  based  on  Dr.  Milano's  resistance  prediction  program 
is  made.  The  resulting  error  in  predicting  resistance  from  errors  made  in  measur¬ 
ing  environmental  and  shipboard  parameters  are  presented. 

The  analysis  indicates  that  velocity,  ice  thickness  and  the  specific 
gravity  of  the  ice  are  the  most  important  parameters  to  be  measured.  Flexural 
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Figure  6. 
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RESISTANCE  (LBS  X  103) 


2 


6 

THICKNESS  (FT.) 
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strength  and  friction  coefficients  must  also  be  measured;  however,  small  errors 
in  these  values  do  not  lead  to  significant  errors  in  the  resistance  predicted. 

Any  errors  made  in  measuring  propulsion  parameters  will  lead  to  a  cor¬ 
responding  error  in  calculating  resistance. 

B.  Introduction 

The  effect  of  various  parameters,  both  environmental  and  machinery,  will 
have  different  impacts  on  computation  of  the  overall  resistance  of  the  POLAR  STAR 
in  ice. 

It  would  be  a  fairly  simple  matter  to  evaluate  the  impact  of  these 
various  parameters  if  there  was  a  verified  algorithm  that  was  available  that 
could  isolate  each  parameter;  unfortunately,  there  is  no  wide  spread  agreement  on 
such  an  algorithm.  Lewis  (1972),  Enkvist  (1972),  Milano  (1972)  and  Vance  (1974) 
have  all  published  information  along  these  lines.  Although  there  is  fair  agreement 
on  some  of  the  parameters  involved,  Milano's  presentation  depends  solely  on  ship 
and  environmental  parameters  and  will  be  used  in  this  study  to  determine  their 
effect  on  the  overall  resistance. 

In  addition  to  the  difficulty  pointed  out  above,  there  is  also  the  prob¬ 
lem  of  the  interaction  of  various  parameters,  i.e.,  the  effect  of  ice  flexural 
strength  increases  with  a  decrease  in  ship  speed  and  ice  thickness.  Thus,  when 
attempting  to  determine  measurement  precision,  these  various  interactions  must  be 
taken  into  consideration. 

This  analysis  is  based  on  the  ship  operating  at  its  design  waterline  of 

28  feet. 

C.  Program  Parameters 

1.  Milano's  Algorithm 

Ship  Parameters  -  Invariant,  determined  from  ship’s  plans 

Environmental  Parameters: 

Thickness  of  ice 
Flexural  strength  of  ice 
Compressive  strength  of  ice 
Density  of  ice 
Density  of  water 
Specific  gravity  of  ice 
Coefficient  of  dynamic  friction 
Specific  gravity  of  water 
Veloci ty 

2.  Vance's  Algorithm 

Ship  Parameters  -  Invariant,  determined  from  ship's  plans 

Environmental  Parameters: 

Thickness  of  ice 

Flexural  strength  of  ice 

Density  of  ice 

Density  of  water 

Coefficient  of  dynamic  friction 

Velocity 
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Eff  ct  of  thickness 


1.  Run  21 

dj:  =  15,000 
V  =  5  FPS  (  3  KTS 


40, 

000 

lbs 

92, 

000 

lbs 

172 

,000 

lbs 

258 

,000 

lbs 

370 

,000 

1  bs 

10 

FPS 

(  6  KTS) 

150,000 

218,000 

310,000 

450,000 

5  FPS  (3  KTS 


50,000  u 
%  INC  TH 


SG 

%  INC  RES 


V  =  10  FPS  (6  KTS! 


^2  2  = 


,000' 

,000. 


R3  0  =  215,000] 
R3’3  =  235,000J 

R4<0  =  315,0001 
R4'4  =  362,000 


16.77 

15.31  MEAN 


9.24  MEAN 


I 


E.  Summary  of  thickness 

1.  At  low  speeds  (5  FPS)  a  5  percent  error  in  ice  thickness  will  lead 
to  a  7  to  8  percent  error  in  resistance. 

2.  At  higher  speeds  (10  FPS)  a  5  percent  error  in  ice  thickness  will 
lead  to  a  5  to  6  percent  error  in  resistance. 

3.  At  low  speeds  (5  FPS)  a  10  percent  error  in  ice  thickness  will 
lead  to  a  15  percent  error  in  resistance. 

4.  At  high  speeds  (10  FPS)  a  10  percent  error  in  ice  thickness  will 
lead  to  a  9  percent  error  in  resistance. 

5.  A  25  percent  increase  in  flexural  strength  has  very  little  effect 
on  the  percent  error  given  above. 

F.  Effect  of  flexural  strength 


1.  Run  32 


o  f  =  12,000 

°c  = 

50,000  m  = 

.2 

SG  = 

.92 

V  =  5  FPS 

V  =  5  FPS 

V  =  10  FPS 

V  = 

10  FPS 

TH  =  2  ft 

TH  =  4  FT 

TH  =  2  FT 

TH 

=  4  ft 

R  =  76,317 

R  =  195,144 

R  =  137,017 

R  = 

249,767 

Run  21 

=  15,000 

°c 

50,000  v  = 

.2 

SG  = 

.92 

V  =  5  FPS 

V  =  5  FPS 

V  =  ;j  FPS 

V  = 

10  FPS 

TH  =  2  ft 

TH  =  4  FT 

TH  -  2  ft 

TH 

=  4  FT 

R  =  92,500 

R  =  258,000 

R  =  150,000 

R  = 

310,000 

Run  30 

of  =  20,000 

°c  = 

50,000  u  = 

.2 

SG  = 

.92 

V  =  5  FPS 

V  =  5  FPS 

V  =  10  FPS 

V  = 

10  FPS 

TH  =  2  FT 

TH  =  4  FT 

TH  =  2  FT 

TH 

=  4  FT 

R  =  127,268 

R  =  386,365 

R  =  173,169 

R  = 

431,200 

Run  27 

,-f  =  30,000 

°c 

50,000  u  = 

.2 

SG  = 

.92 

V  =  5  FPS 

V  =  5  FPS 

V  =  10  FPS 

V  = 

10  FPS 

TH  =  2  FT 

TH  =  4  FT 

TH  =  2  FT 

TH 

=  4  FT 

R  =  214,812 

R  =  731,220 

R  =  234,237 

R  = 

758,726 
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5.  RUN  36 


0 


of  =  14,300 

°c 

50,000  u  = 

.2  SG  =  . 

92 

V  =  5  FPS 

V  =  5  FPS 

V  =  10  FPS 

V  =  10  FPS 

TH  =  2  FT 

TH  =  4  FT 

TH  =  2  FT 

TH  =  4  FT 

R  =  89,165 

R  =  242,246  R  =  146,192 

R  =  294,426 

6.  RUN  37 

o  ^  —  16,000 

°c 

50,000  m  = 

.2  SG  =  . 

92 

V  =  5  FPS 

V  =  5  FPS 

V  =  10  FPS 

V  =  10  FPS 

TH  =  2  FT 

TH  =  4  FT 

TH  =  2  FT 

TH  =  4  FT 

R  =  99,607 

R  =  381,174  R  =  153,616 

R  =  331,355 

Calculation  Data 

%  INC  g f 

%  IN  R 

SPEED* 

THICKNESS* 

5%  (14.3K  to  15K) 

3.6 

L 

L 

6.1 

L 

H 

2.5 

H 

L 

5.0 

H 

H 

10%  (14.3K  to  16K) 

10.0 

L 

L 

13.8 

L 

H 

4.8 

H 

L 

11.14 

H 

H 

20%  (12K  to  15K) 

17.5 

L 

L 

24.4 

L 

H 

8.7 

H 

L 

19.4 

H 

H 

25%  ( 15K  to  20K) 

27.3 

L 

L 

33.2 

L 

H 

13.4 

H 

L 

28.1 

H 

H 

*L  in  speed  column 

refers  to  the 

low  speed  range. 

H  refers  to  the 

higher 

speeds  used  in  the  sample  calculations.  The  same 

is  true  for  the 

thick- 

ness  used.  The  exact  number  can  be  found  by  referring  to  the  above 
sample  calculations. _ 


G.  Summary  of  flexural  strength 

The  effect  of  flexural  strength  is  related  to  the  speed  of  the  vessel  and 
the  thickness  of  the  ice.  As  shown  in  the  accompanying  example  calculations,  the 
greatest  effect  occurs  at  low  speed,  high  thickness.  In  this  case  a  5  percent 
error  in  measuring  flexural  strength  will  lead  to  a  6.1  percent  error  in  resis¬ 
tance.  At  high  speeds  with  low  thicknesses  the  results  show  an  opposite  trend, 
i.e.,  5  percent  error  in  measuring  flexural  strength  will  lead  to  a  2.5  percent 
error  in  resistance. 

It  should  be  noted  that  the  results  presented  here  are  from  Milano's 
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(1972)  program  and  show  a  slightly  stronger  dependence  of  resistance  and  flexural 
strength  than  shown  by  Lewis  (1972)  and  Vance  (1974).  This  fact  should  be  kept 
in  mind  when  utilizing  the  data  presented  herein. 

H.  Effect  of  ice  specific  gravity 

1.  RUN  28 


af  -- 

:  15,000 

a  =  50 
c 

,000 

u  - 

.2 

SG 

V  = 

5  EPS 

V  = 

5  FPS 

V  = 

10 

FPS 

V  = 

10 

FPS 

TH  = 

:  2  FT 

TH 

=  4  FT 

TH  = 

=  2 

FT 

TH 

-  4 

FT 

R  = 

106,971 

R  = 

304,183 

R  = 

162 

,786 

R  = 

355 

,560 

RUN 

21 

cf  = 

15,000 

o„  =  50 
c 

,000 

u  = 

.2 

SG 

V  = 

5  FPS 

V  = 

5  FPS 

V  = 

10 

FPS 

V  = 

10 

FPS 

TH  = 

2  FT 

TH 

=  4  FT 

TH  = 

=  2 

FT 

TH 

=  4 

FT 

R  = 

92,500 

R  = 

258,000 

R  = 

150 

,000 

R  = 

310 

,000 

RUN 

33 

°f  = 

‘  15,000 

o  =  50 
c 

,000 

u  = 

.2 

SG 

V  = 

5  FPS 

V  = 

5  FPS 

V  = 

10 

FPS 

V  = 

10 

FPS 

TH  = 

2  FT 

TH 

=  4  FT 

TH  = 

=  2 

FT 

TH 

=  4 

FT 

R  = 

86,900 

R  = 

235,722 

R  = 

142 

,715 

R  = 

287 

,099 

Calculation  Data 


%  INC  SG 

1  DEC  R 

SPEED 

THICKNESS 

10%  (.82  to  .92) 

15.6 

L 

L 

17.9 

L 

H 

8.5 

H 

L 

14.7 

H 

H 

5%  (.96  to  .966) 

6.4 

L 

L 

9.4 

L 

H 

5.1 

H 

L 

8.0 

H 

H 

15S  (.82  to  .966) 

23.1 

L 

L 

29.0 

L 

H 

14.1 

H 

L 

23.8 

H 

H 

I.  Summary  of  specific  gravity 

Although  the  specific  gravity  of  the  ice  falls  within  small  limits,  i.e., 
.89  to  .94  with  a  mean  of  .92,  errors  made  in  measuring  the  S.G.  will  have  a  large 
effect  on  the  calculation  of  the  resistance  of  the  vessel.  This  is  particularly 
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true  at  low  speeds  and  high  thicknesses,  which  one  would  suspect,  for  it  is  in 

this  operational  mode  that  the  energy  utilized  to 

submerge  the  ice  is  the  greatest 

percentage  of  the  overall  resistance.  However,  one  should 

be  aware,  the  thicker 

the  ice  the  greater  the  impact  of 

the  specific  gravity  of 

the  ice. 

The  sample  calculations 

using  Milano's 

program  indicate  a  5  percent 

error  in  S.6.  will  lead  to  a  9.4 

percent  error  in 

resistance  at  the  worst  condition 

A  15  percent  error  in  S.G.  will  lead  to  a  20  percent  error 

in  resistance  at  the 

worst  condition. 

J.  Effect  of  velocity 

1.  Run  21 

af  =  15,000 

oc  =  50,000 

u  = 

.2  SG  =  .92 

TH  =  1  FT 

%  INC  VEL 

;  INC  R 

5-7  FPS 

30 

36 

7.5-10  FPS 

25 

40 

9.5-10  FPS 

10 

18 

9.5-10  FPS 

5 

6.4 

TH  =  2  FT 

5  -  7 

30 

16.3 

7.5  -  10 

25 

23.3 

9  -  10 

10 

12.0 

9.5  -  10 

5 

6.7 

TH  =  4  FT 

i  5  -  7 

30 

6.5 

7.5  -  10 

25 

9.0 

9  -  10 

10 

2.9 

9.5  -  10 

5 

1.6 

TH  =  6  FT 

5  -  7 

30 

3.1 

7.5  -  10 

25 

5.3 

9  -  10 

10 

2.1 

9.5  -  10 

5 

1.3 

2.  RUN  30 

-f  =  20,000 

■jc  =  50,000 

u  - 

.2  S3  =  .92 

TH  =  1  FT 

%  INC  VEL 

INC  R 

TH  =  2  FT 

%  INC  VEL 

%  INC  R 

5  -  7 

30 

9.4 

7.5  -  10 

25 

16.6 

9  -  10 

10 

7.9 

9.5  -  10 

5 

3.9 

TH  =  4  FT 

5  -  7 

30 

4.4 

7.5  -  10 

25 

4.9 

9  -  10 

10 

.9 

9.5  -  10 

5 

.5 

TH  =  6  FT 

5  -  7 

30 

10 

7.5  -  10 

25 

3.4 

9  -  10 

10 

1.5 

9.5  -  10 

5 

.7 

K.  Summary  of  velocity  effect 

The  effect  of  velocity  is  more  pronounced  at  low  thickness  and  in  we.'  ^r 
ice  as  one  would  suspect.  The  sample  calculations  indicate  a  10  percent  error  in 
velocity  measurements  will  lead  to  an  18  percent  error  in  resistance  in  one  foot 
of  ice.  In  2  feet  of  ice  a  10  percent  error  in  velocity  will  lead  to  a  12  percent 
error  in  resistance.  In  6  feet  of  ice  when  other  resistance  components  predomi¬ 
nate,  a  10  percent  error  in  velocity  measurement  will  only  cause  a  2  percent  error 
in  resistance. 

The  interaction  of  the  various  components  can  be  seen  from  the  decrease 
in  velocity  effect  where  the  flexural  strength  is  increased.  The  25  percent  in¬ 
crease  in  strength  causes  a  decrease  in  the  effect  of  velocity  on  the  resistance. 

The  importance  of  accurately  measuring  velocity  is  not  only  emphasized 
by  the  errors  pointed  out  herein,  but  in  addition,  the  velocity  of  the  ship  is 
utilized  directly  when  making  power  calculations  from  machinery  readouts. 

L.  Effect  of  compressive  strength  (ac) 

1.  Run  35 

af  =  15,000  oc  =  45,000  u  =  .2  SG  =  .92 

V  =  5  FPS  V  =  5  FPS  V  =  10  FPS  V  =  10  FPS 

TH  =  2  FT  TH  =  4  FT  TH  =  2  FT  TH  =  4  FT 

R  =  94,579  R  =  261,071  R  =  105,394  R  =  312,448 

2.  Run  21 

af  =  15,000  oc  =  50,000 

V  =  5  FPS  V  =  5  FPS  V  =  10  FPS  V  =  10  FPS 

TH  =  2  FT  TH  =  4  FT  TH  =  2  FT  TH  =  4  FT 

R  =  93,000  R  =  258,000  R  =  150,000  R  =  310,000 
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3.  RUN  34 


o  ^  =  15,000 


ac  =  52,500 


V  =  5  FPS 

V  = 

5  FPS 

V  =  10 

FPS  V  =  10  FPS 

TH  =  2  FT 

TH 

=  4  FT 

TH  =  2 

FT  TH  =  4  FT 

R  =  92,829 

R  = 

256,421 

R  =  148 

,644  R  -  307,798 

4.  RUN  17 

=  15,000 

ac  =  100,000 

V  =  5  FPS 

V  = 

5  FPS 

V  =  10 

FPS  V  =  10  FPS 

TH  =  2  FT 

TH 

=  4  FT 

TH  =  4 

FT  TH  =  4  FT 

R  =  86,813 

R  = 

240,437 

R  =  142 

,628  R  =  291,814 

Calculation  Data 

RUN  S 

IN  C 

ac 

%  INC  R 

SPEED 

THICKNESS 

(34-21) 

5% 

.2 

L 

L 

.6 

L 

H 

.9 

H 

L 

.7 

H 

H 

(21-35) 

10 % 

1.7 

L 

L 

1.2 

L 

H 

.3 

H 

L 

.8 

H 

H 

,  (34-35) 

15% 

1.8 

L 

L 

1.8 

L 

H 

1.2 

H 

L 

1.5 

H 

H 

(21-17) 

50% 

6.6 

L 

L 

6.8 

L 

H 

4.9 

H 

L 

5.9 

H 

H 

M.  Summary  of  compressive  strength  effect 

Small  changes  in  the  compressive  strength  of  the  ice  has  very  little  effect 
on  the  resistance  predicted.  It  is  not  until  we  encounter  errors  in  the  range  of 
50  to  60  percent  do  we  see  a  5  to  6  percent  change  in  resistance. 

For  all  practical  purposes,  acceptable  values  of  compressive  strength  of 
ice  presented  in  the  literature  can  be  used  without  any  significant  effect  on  the 
resistance  prediction. 

N.  Effect  of  friction  (u) 
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msm 


V  =  5  FPS 

TH  =  2  FT 

V  =  5  FPS  V  = 

TH  =  4  FT  TH  = 

10  FPS 
=  2  FT 

V  =  10  FPS 
TH  =  4  FT 

R  =  78,533 

R  =  206,913  R  = 

140,948 

R  =  262,825 

RUN  21 

=  15,000 

oc  =  50,000 

u  =  .2 

SG  = 

V  =  5  FPS 

TH  =  2  FT 

R  =  92,500 

V  =  5  FPS  V  = 

TH  =  4  FT  TH= 

R  =  258,000  R  = 

10  FPS 

2  FT 
150,000 

V  =  10  FPS 
TH  =  4  FT 

R  =  310,000 

RUN  40 

=  15,000 

ac  =  50,000 

u  =  .1 

SG  = 

V  =  5  FPS 

TH  =  2  FT 

R  =  100,786 

V  =  5  FPS  V  = 
TH  =  4  FT  TH  = 
R  =  283,325  R  = 

10  FPS 
;  2  FT 
153,301 

V  =  10  FPS 

V  =  4  FT 

R  =  332,435 

RUN  41 

o.p  =  15,000 

oc  =  50,000 

u  =  .235  SG  =  . 

V  =  5  FPS 

TH  =  2  FT 

R  =  98,561 

V  =  5  FPS  V  = 
TH  =  4  FT  TH  = 
R  =  275,684  R  = 

10  FPS 

2  FT 
152,066 

V  =  10  FPS 
TH  =  4  FT 

R  =  325,474 

RUN  39  (.1)  - 
RUN  21  (.2)  - 

21  (.2)  50%  RUN  21  (.2)  - 
40  (.25)  20%  RUN  41  (.235) 

41  (.235) 

-  40  (.25) 

Calculation  Data 

%  INC  u 

%  INC  R 

SPEED 

THICKNESS 

6%  (41-40) 

2.2 

L 

L 

2.7 

L 

H 

.8 

H 

L 

2.1 

H 

H 

15%  (21-41) 

6.1 

L 

L 

6.4 

L 

H 

1.4 

H 

L 

4.8 

H 

H 

20%  (21-40) 

8.2 

L 

L 

8.9 

L 

H 

2.1 

H 

L 

6.7 

H 

H 

5 


4 


Calculation  Data  (Continued) 
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%  INC  n 

%  INC  R 

SPEED 

THICKNESS 

50%  (39-21) 

15.1 

L 

L 

19.8 

L 

H 

6.0 

H 

L 

15.2 

H 

H 

O.  Summary  of  friction  effect 

The  effect  of  the  coefficient  of  friction  (defined  as  the  normal  force 
divided  by  the  tangetral  force,  for  the  lack  of  any  better  understanding  of  the 
phenomenon)  is  difficult  to  appreciate.  If  one  is  dealing  with  dynamic  friction, 
the  coefficient  can  vary  from  .025  to  as  high  as  .5.  Most  of  the  published  data 
confines  the  range  from  .1  to  .3. 

Milano's  program  indicates  that  for  small  changes  in  the  friction  co¬ 
efficient,  i.e.,  6  percent  (.235  to  .25)  there  is  only  a  small  effect  on  the  over¬ 
all  resistance,  i.e.,  3  percent.  However  a  50  percent  increase,  .1  to  .2,  could 
cause  a  70  percent  error  in  the  resistance  estimate. 

Recent  model  tests  (Lecourt,  1975)  have  indicated  that  an  increase 
in  the  coefficient  from  .037  to  .28  (7.5  fold  increase)  has  lead  to  a  40  percent 
increase  in  resistance.  Full  scale  tests  run  by  the  Wartrila  Shipyard  in  Finland 
(Makimer,  1975)  have  indicated  a  30  percent  decrease  in  resistance  of  the  ice¬ 
breaking  vessel  SILMA  after  being  coated  with  a  low  friction  coating.  Unfortunately, 
no  friction  measurements  were  taken  on  the  vessel. 

In  summary,  it  can  be  said  that  small  errors  in  measuring  the  friction 
coefficient,  i.e.,  >10  percent,  will  have  only  a  minor  effect  on  the  resistance 
prediction,  however,  one  must  know  or  have  some  good  idea  of  what  the  friction 
coefficient  is,  i.e.,  .10  or  .20. 

P.  Effect  of  propulsion  equipment  measurements 

From  the  equations  in  the  next  section,  it  can  be  seen  that  the  relation¬ 
ship  between  the  measured  parameters  and  the  computed  resistance  in  a  linear  one 
and  any  error  made  in  measuring  the  parameter,  i.e.,  torque,  RPM,  velocity  or 
propulsive  coefficient  (which  is  estimated  from  the  model  tests)  will  result  in  an 
error  of  equal  magnitude  in  computing  the  resistance.  Several  examples  of  typical 
values  are  presented. 

Q.  Equations 

1.  Effect  of  thrust 

RESISTANCE  =  THRUST^  SHAFT)  x  P.C. 

@  13  FPS  4'  THICKNESS  R  =  338,633 


THRUST  =  —5—  =  =  564,388 

P.C.  .6 


« 


1 


52 


'SHAFT 


=  564,388 


%  INC  T 


SHAFT 


%  INC  R 


P.C.  =  .6 
RESIS  =  338,632 


(T  =564388-T= 536000) 


'SHAFT 


=  536,000 


P.C.  =  .6 
RESIS  =  321,600 


(T=564388-T= 507000) 


'SHAFT 


=  507,000 


P.C.  =  .6 
RESIS  =  304,200 


'SHAFT 


=  564,388 


P.C.  - 


(P.C.=.6=P.C.+.57) 


RESIS  =  338,632 


'SHAFT 


P.C.  = 


=  564,388 


(P.C.=.6-P.C.=. 54) 


RESIS  =  321,701 


‘SHAFT 
P.C.  = 


=564,388 


RESIS  =  304,769 


2.  Effect  of  torque 
_  R  n  Q  N 


EHP  =  SHP  x  P.C. 


RESISTANCE  = 


EHP  x  550 


RESISTANCE  =  2  n  QN  -P-C-: 


53 


9 


FACTOR/ERROR 


FACTOR/ ERROR 


1 


J 


s 


I 


I 


I 


« 


Vel .  - 

5-7  FPS 
7.5  -  10  FPS 
9  -  10  FPS 
9. *  -  10  FPS 


P.C.  = 

.6 

N  = 

2.92 

RPS 

5% 

.57 

5% 

2.77 

RPS 

10% 

.54 

10% 

2.63 

RPS 

RUN  21  13  FPS 


ICE:  EHP 


RES  x  VEL  338633  x  .3 
550  “  550 


=  8004 


SHP 


EHP 

P.C. 


=  -8°°4  =  13340 


n  _  550  x  SHP 
g  "  2  ::  N 


=  399904 


Q  = 

399904 

%  INC  Q 

N  = 

2.92  RPS 

5 

V  = 

5  FPS 

(Q- 399904-Q=380000) 

P.C.  = 

.6 

R  = 

880428 

10 

Q  = 

380000 

(Q=399904-Q=360000) 

N  = 

2.92  RPS 

V  = 

5  FPS 

P.C.  = 

.6 

R  = 

836608 

Q  = 

360000 

N  = 

2.92  RPS 

V  = 

5  FPS 

P.C.  = 

.6 

R  = 

792576 

Q  = 

399904 

%  INC  N 

N  = 

2.92  RPS 

5 

V  = 

5  FPS 

(N=2.92  -  N=2.77) 

~o 

o 

II 

.6 

R  = 

880438 

10 

Q  = 

399904 

( N=2 . 92  -  N=2. 63) 

N  = 

2.77  RPS 

V  = 

5  FPS 

P.C.  = 

.6 

R  = 

835210 

Q  = 

399904 

N  = 

2.63  RPS 

V  = 

5  FPS 

P.C.  = 

.6 

R  = 

792997 

Q  = 

399904 

%  INC  P.C. 

N  = 

2.92 

5 

V  = 

5  FPS 

(P.C. =.6  -  P.C. =.57) 

P.C.  = 

.6 

R  = 

880439 

10 

Q  = 

399904 

(P.C. =.6  -  P.C. =.54) 

N  = 

2.92 

V  = 

5  FPS 

54 


%  INC  R 
4.9 


9.9 


%  INC  R 
5.1 


9.9 


%  INC  R 
5 


10 


■ 

i, 


i 


P.C.  =  .57 

R  =  836417 
Q  =  399904 
N  =  2.92 
V  =  5  FPS 
P.C.  =  .54 

R  =  792395 

n  _  2nQN  P.C. 
K  "  V 


Q  =  399904 

N  =  2.92  P.C 

.  =  .6 

! 

VEL 

%  INC  VEL 

RESISTANCE 

%  INC  R 

5 

7 

30 

880439 

628885 

28.5 

7.5 

10 

25 

586959 

440219 

25 

< 

9 

10 

10 

489133 

440219 

10 

9.5 

10 

5 

463389 

440219 

5 

R.  Electrical 

power  measurement 

By  noting  voltage  and  amperage 

to  the  motor, 

a  rough  estimate 

power  being  used  can 

be  obtained. 

P  =  E  I  watts 
SHP  =  watts/746 

RESISTANCE  =  -HP  *  550 

El  x  550  P.C. 


746  V 

RESISTANCE  =  ~37?,  EI  P'C' 


Thus,  as  in  the  previous  equations,  due  to  the  linear  relationship 
between  power  and  resistance,  an  error  in  measuring  voltage  or  amperage  would 
lead  to  an  equal  percentage  error  in  the  resistance  computed. 

S.  Summary  of  propulsion  equipment 

It  is  evident  from  the  equations  presented  here  and  in  other  sections 
of  this  report  that  the  relationship  between  propulsion  equipment  measurements, 
i.e.,  torque,  thrust,  volts,  amps,  RPM  and  the  prediction,  or  in  this  case,  the 
calculation  of  resistance,  is  a  direct  one.  That  is,  any  error  made  in  measuring 
these  parameters  leads  to  an  error  of  equal  magnitude  in  the  resistance  calculated. 


T.  Effect  of  draft 

1.  RUN  45  of  =  15,000 


c  =  50,000 


=  .2 


SG  -  .92  DR  28 


55 


2.  RUN  44 


Of  =  15,000 


oc  =  50,000 


p  =  .2 


SG  =  .92  DR  30 


THICKNESS 

=  2 

FT 

VEL 

RESISTANCE  (45) 

RESISTANCE  (44) 

%  DIFF 

5 

95,440 

92,088 

3.5 

10 

149,259 

144,460 

3.2 

15 

261,422 

253,155 

3.1 

20 

425,813 

411,587 

3.3 

THICKNESS 

=  4 

FT 

5 

266,586 

251,771 

5.5 

10 

313,702 

293,055 

6.6 

15 

369,833 

345,800 

6.4 

20 

486,721 

454,911 

6.5 

THICKNESS 

=  6 

FT 

5 

585,930 

558,699 

4.6 

10 

642,720 

611,227 

4.8 

15 

698,052 

656,581 

5.9 

20 

771,628 

724,291 

6.1 

U.  Summary  of  draft  effect 

The  change  in  draft  of  the  POLAR  STAR  from  28  feet  to  30  feet  changes 
many  of  the  parameters  that  enter  the  icebreaking  prediction  equation,  i.e.,  dis¬ 
placement  changes  from  10,860  tons  to  12,200  tons,  the  bow  angle  change  from  15° 
to  14.5',  etc.  Therefore,  a  whole  new  data  set  was  used  in  Dr.  Milano's  program 
for  a  draft  of  30  feet  (see  Appendix  B).  The  results  of  the  comparison  indicate 
that  when  the  draft  is  increased  to  30  feet,  the  overall  ice  resistance  decreased 
from  3  to  6  percent  with  the  larger  percentages  at  the  higher  ice  thicknesses. 

See  Figure  8. 

In  addition  to  the  icebreaking  capability,  there  are  two  other  aspects 
of  increased  draft  that  must  be  kept  in  mind.  Once  is  the  effect  on  the  extractor 
difficulty  encountered;  this  effect  will  be  covered  in  the  Ramming  and  Extractor 
Prediction  aspect.  The  record  aspect  is  the  decreased  clearance  (overhang)  caused 
by  increased  draft.  The  POLAR  STAR  was  designed  with  a  vertical  clearance  of  5 
feet  between  the  28  foot  waterline  and  the  start  of  the  curvature  of  the  stern  to 
the  vertical.  The  increased  draft  would  decrease  this  clearance  to  approximately 
3  feet. 

VI.  RAMMING  AND  EXTRACTION  PREDICTIONS  FOR  THE  POLAR  STAR  BY  CDR  G.  P.  VANCE 

A.  Abstract 

This  report  contains  predictions  for  the  POLAR  STAR  in  the  ramming  mode. 
The  calculations  indicate  that  the  POLAR  STAR  would  be  capable  of  breaking  almost 
21  feet  of  homogenous  ice  in  ramming  (see  Appendix  C). 

Utilizing  the  theory  available  the  extraction  difficulty  is  also  pre¬ 
dicted. 

B.  Introduction 

The  prediction  of  a  vessel's  capability  to  break  ice  in  the  ramming  mode, 
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RESISTANCE  (LBS  X  10J) 


i,e.,  backing  up  and  ramming  into  the  ice  and  coming  to  a  complete  stop,  has  been 
attempted  by  White  (1969)  and  modified  by  Lewis  (1970).  The  method  assumes  a  rigid 
flat  sheet  of  homogenous  ice  with  a  specified  tensile  strength.  The  vessel  rides 
up  on  the  ice  sheet  and  exerts  a  downward  force  on  the  ice  until  failure  occurs. 

The  downward  force  generated  is  a  function  of  the  surge,  heave  and  pitch  motions 
of  the  ship  and  the  impact  velocity  and  displacement  of  the  ship  (see  Fig.  9). 

White  assumed  a  realtively  simple  relationship  between  the  downward  force 
exerted  to  break  the  ice  and  Lewis  introduced  some  field  data  to  refine  the  force 
equation. 

The  extraction  difficulty  is  a  more  elusive  number  to  pin  down.  White 
assumed  that  the  extraction  difficulty  was  a  function  of  the  downward  force,  the 
coefficient  of  static  friction  and  the  bow  angles.  He  then  formed  a  ratio  of  the 
thrust  needed  to  extract  the  vessel  from  the  ice  and  the  bollard  foreward  thrust 
available  from  the  vessel.  This  "Difficulty  of  Extraction  Ratio"  was  not  an 
absolute  number,  i.e.,  it  was  a  number  one  could  compare  to  other  vessels  such  as 
the  WIND  class  and  the  GLACIER. 

This  report  presents  calculations  made  for  the  POLAR  STAR. 

Where:  B  =  maximum  beam  in  feet 

d  =  draft  of  the  vessel  in  feet 
=  water  plane  area  coefficient 
Cb  =  block  coefficient 
SA  =  spread  angle  complement  in  radians 
f|<  =  coefficient  of  kinetic  friction  (assumed  to  be  .2  for 
uni formi ty ) 

BA  =  bow  angle  in  radians 

Figure  9  is  a  plot  of  the  predicted  ice  thickness  that  can  be  broken  by 
the  POLAR  STAR.  The  sample  calculations  are  shown  in  the  appendices.  It  should 
be  noted  the  plot  reflects  a  28  foot  draft  and  f  =  235  PSI. 

Lewis  (1970)  has  indicated  that  k  =  .45  used  by  White  may  have  been  too 
low  and  recommends  a  figure  of  k  =  1.02  from  full  scale  tests,  he  also  indicates 
that  a  f  =  235  PSI  may  be  too  high  and  a  more  realistic  figure  would  be  in  the 
area  of  100  to  140  PSI . 

C.  Ramming  capability 

The  equations  used  in  predicting  the  thickness  of  ice  by  White  are: 

F  =  k  t^ 
rBZR  K  L 

substituting  value  of  .45  for  k  and  235  PSI  for  the  tensile  (flexural)  strength  of 
ice: 


tr  FBZR 
100 

where  "tr"  stands  for  thickness  of  standard  ice  broken  due  to  ramming,  in  inches, 
and  "FgR7"  is  the  relatively  sustained  downward  force  against  the  ice  in  pounds. 
This  last  parameter  "Fgp/1  further  depends  on  three  factors:  impact  velocity  "v" 
in  FPS,  ship  displacement  "A"  in  pounds,  and  the  White  Ratio  "WR".  The  equation 
for  calculating  "Fg^"  is  as  follows: 

FBZR  =  6-64  \l.  (WRXa)0-845 
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PREDICTION  OF  THICKNESS  OF  STANDARD  SEA  ICE  BROKEN 
AS  A  RESULT  OF  RAMMING  AS  A  FUNCTION  OF  IMPACT  VELOCITY 


LENIN 

POLAR  STAR 
ST.  LAURENT 

GLACIER 


WIND 

(WITHOUT  BOSSING) 


The  White  Ratio  is  a  coefficient  of  performance  that  was  developed 
by  White,  in  his  study  of  the  mechanics  of  icebreaking.  The  equation  used  to 
obtain  this  numbers  is: 

WR  =  1. 000234(10. 72+B/d)(0.1833+CW)(1.652-Cb)  x 
(6. 14-SA2) (0. 725-fk) ( 1. 718-BA) 

The  recommended  equations  would  be: 

fbrz  =  1-02  °ft2  or 

tr  =  FBZR 
106.25 

i.e.,  of  104.16  PSI  =  15000  PSF  =  7.32  kg/cm2 

D.  Extraction  difficulty 

White  (1969)  has  proposed  an  index  of  extraction  difficulty  that  reflects 
a  computed  extraction  thrust  based  on  the  hull  form  and  coefficient  of  static 
friction  divided  by  the  maximum  bollard  thrust  available.  This  ratio  is  more  of 
a  comparison  to  other  vessels  rather  than  an  exact  indication  of  the  thrust  avail¬ 
able  or  necessary  to  extract  the  vessel  after  ramming.  Figure  10  shows  the  com¬ 
puted  index  of  extraction  for  several  icebreakers  and  that  of  the  POLAR  STAR.  The 
POLAR  STAR  curve  is  presented  for  a  bow  angle  of  15°  and  a  bow  angle  of  30  .  The 
equations  are  very  sensitive- to  the  bow  angle  and  this  fact  is  reflected  in  the 
careful  design  of  the  POLAR  STAR  bow.  The  POLAR  STAR  has  a  bow  angle  of  15;  at 
the  28  foot  water  line  this  angle  gradually  increases  to  30°  at  the  18  foot  water 
line  where  it  is  carried  to  a  6  foot  vertical  step.  Thus  when  computing  the  ex¬ 
traction  difficulty  one  must  determine  how  far  the  ship  will  penetrate  the  ice  to 
determine  its  extraction  difficulty.  In  most  instances  of  ramming  the  momentum  of 
the  vessel  will  be  sufficient  to  carry  it  to  a  trim  of  1°  or  2°  which  will  insure 
a  bow  angle  of  greater  than  15°  and  close  to  30°.  Thus  the  POLAR  STAR  should  not 
encounter  any  unusual  difficulty  in  extracting  itself  from  a  ram  (Fig.  10). 


Figure  10 
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APPENDIX  A 

Curves  of  Resistance  versus  Velocity 
and  Thrust  versus  Velocity 


RESISTANCE  (LBS  X  103) 


2  4  6  8 

THICKNESS  (FT.) 

Figure  15 
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APPENDIX  B 

POLAR  STAR  30  ft  Draft  Particul 


POLAR  STAR  30'  particulars 

LCF  =  2.4  ft  fsd 

LWL  =  363.75  ft 

LCG  =  3.1  ft  fwd 

d  =  30  ft 

K  6  =  28  ft 

BWL  =  79.33  ft 

=  12,200  tons  -  27,328,000  lbs 
BA  =  14. 5-  =  25307  RAD 
SA  =  .42  RAD 

Fs  =  -6 
Cw  =  .762 

CM  =  .85 

Cp  =  .585 

KMl  =  379.5 

GMl  =  351.5 

Cb  -  .  5 

Water  line  entrace  angle  =  59.5° 

WR30  =  (.000234)  (10. 72+  ^^)  ( .  1833+.  762)  ( 1. 652- .  5) 

(6.14  -  (.42)2)  (. 725-. 2) ( 1. 7180. 25307) 

=  (.000234) (13. 364) (.9453) (1.152) (5. 9636) (.525) (1.457) 
=  .0155346 

FBZR  =  6.64  V  (.0155346  x  2. 7328xl07) *845 
=  6.64  V  ( 5. 69617xl04) 

=  3. 78226xl0F  (V) 
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APPENDIX  C 


SAMPLE  CALCULATIONS 


POLAR  STAR 


Sample  calculations  of  ramming  capability 
(28  foot  draft) 


100 

2.  Fbzr  =  6.64  v  (WRxa)0-845 

3.  WR  =  0. 000234(10. 72+B/d)(0.1833+Cu)(1.652-Ch)x 

(6- 14-SA2)(0.725-fk)(1.718-BA)  D 

POLAR  STAR 

B  =78.0  feet 
d  =  28.0  feet 
Cw  =  .74 
Cb  =  .49 
SA  =  .420 
fk  =  .2 
BA  =  .261 

=  24,333,120  pounds 

First,  the  White  Ratio,  WR  was  computed: 

WR  =  0.000234(10.72  +78- °/28.0) (0. 1833  +  . 74) ( 1.652-. 49)x 
6.14-(.420)2  (0.725- . 2) ( 1. 718-.261) 

WR  =  .015467 


Then  FRZR  @  speeds  of  5  FPS ,  10  FPS,  15  FPS,  and  20  FPS  were  calculated  in 
the  respective  order: 


@5FPS : 

@10FPS : 

015FPS: 

@20FPS : 

Finally  tr 
05FPS : 

010FPS 


fBZR  =  6-64(5) ( . 01546x24, 333, 120)0'845 
=  1,707,486.35  pounds 

FR7R  =  6. 64(10)(. 01546x24, 333, 120)0*845 
=  3,414,972.71  pounds 

Fbzr  =  6. 64( 15) ( . 01546x24, 333, 120)0*845 
=  5,122,459.06  pounds 

Fbzr  =  6. 64(20) ( .  01546x24, 333, 120)°‘845 


=  6,829,945.42  pounds 


was  computed: 

*  _  1,707,486.35 

V  "  ioo 

*  3,414,972.71  = 

V  "  ioo 

=  5,122,459.06  = 


•  =  130"  or  10.8' 
184"  or  15.4' 

226"  or  18.9' 


0 1 5  FPS  t 


I 


@20FPS  tr  =  6-8iqq9-4-5-  =  261"  or  21.8' 
The  following  is  a  table  of  values  computed 


V 

fbrz 

tr  in. 

tr  ft 

5 

1,707,486.35  pounds 

130 

10.8 

10 

3,414,972.71  pounds 

184 

15.4 

15 

5,122,459.06  pounds 

226 

18.9 

20 

6,829,945.42  pounds 

261 

21.8 

The  values  of  tr  in  inches  were  then  plotted  against  a  scale  of  impact 
velocity  in  FPS  using  speeds  of  5,  10,  15,  and  20  see  Figure  9. 

POLAR  STAR 

INDEX  OF  EXTRACTION  DIFFICULTY  (EQNS) 

CALCULATION  OF  Et 
Extraction  difficulty  index  =  Et 

t: 


L  Ft  =  fe 


;os --sine 


2.  a7  =  (cosSA)  cos(BA+  )  +  f.  sin(BA+o) 


3.  a7  =  (cosSA)  sin(BA+  )  +  f  cos(BA+o) 

4.  =  A4C3  "  A3C4 


A4C3  "  ^3C4 
^3^4  “  A4O3 


5.  A3  =(LBP)(B)(CW)(C4.2)  (Constant) 

6-  A4  =  FBZR 

tan  BA 

7.  B3  =  A3  (LCG  -  LCF)  (Constant) 

S.  B4  =  Fbzr  (d-KG)  +  LBP  -  LCG  _  (a_fbzr)  (GMj) 


tanBA 


9-  C3  =FBZR 


!0.  C4  -Fbzr  LBP  _  LCG 


« 


74 


POLAR  STAR  28'  particulars 


LCF  = 

1.9  ft 

BA  =  .261 

LBP  = 

352.0  ft 

SA  =  .420 

LCG  = 

3.3  ft 

fs  =  .2  x  .6 

d  = 

28.0  ft 

CW  =  .74 

KG 

28.0  ft 

CB  =  -49 

B  = 

78.0  ft 

KM(_=  378.1  ft 

A  - 

1086 3 ( tons ) x2240( 1 bs )  = 

6Ml=350. 1 

24333120  lbs 


POLAR  STAR 

Index  of  extraction  difficulty  (calculations) 

Calculation  of 

•  @  5FPS ,  FgzR  =  1,707,486.354 
A3  =  ( 352. 0) ( 78.0)( . 74) (64.2)  =  1,304,379.65 

A4  =  ~ =  6392862.96 

B3  =  (1304379. 65) ( 3. 3-1. 9)  =  1,826,131.51 

B.  =  -1707486.354  (28-28)  +  -  3.3  -  (24,333,120-1707,486.354) 

Lt26T)  (35°-1) 

=  1707486.354  ( 646 . 8 )- 792 1234339 
B4  =  1104430312  x  -9,025,636,513 

C3  =  1707486.354 

C,  =  (1707486.354)  352-3.3  =  294,882,893.3 
4  2 


A3B4-A4B3 

=  (6,392.863)  ( 1 ,707 ,486 . 354)- ( 1 ,304 ,379.65) (294 ,882 ,893) 
(  1,304, 379. 65) (9, 025, 636. 513) -(6, 392, 863) (  1  ,826, 131) 

=  (1.0916xl013)  -  3.846xl014 
(-1. 17728xl016)-l. 1674x1013 

=  -3.73681xl014  =  .03172  Radians 
-  1.178x  1016 

a-j  =  cosSA  cos(BA+c)  +  f$  sin  (BA+n) 
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=  cos  .42  cos  (.262+. 03172)+. 6  sin  (. 262+. 03172) 
=  ( . 913) ( . 957)  +  .6  (.289) 


=  .874  +  .1734 
=  1.0474 

b i  =  (-cosSA)  sin  (BA+o)  f$  cos  (BA+o) 
=  (-.913)  (.2895)  +  .6  (.957) 


-.2643  +  .5743  - 

1.0474 


=  3.378 


.31 


.31 


FBRZ 

a7 

r -  COS') 


si  no 


=  1. 707486xl06 

3. 378(cos  .03172)  -  sin  .03172 


=  1.707486x10° 
3. 378( . 99949 


1. 707486xl06 
3.3763  -  .031 


=  5.10522xl05 

T1D  =  973  700  lbs  (LASKY,  NSRDC  Rpt,  1971) 
Id 

=  5.10522xl05  =.52431 

IB6  9.737x105 

•  0  10FPS ,  FgRZ  =  3,414973  lbs 
A3  =  1304378 


B3  =  1826132 

C3  =  Fbzr  =  3.414973 

A4  =  [BZR  =  3414973  =  1.27857xl07 

tan  BA  tan. 261 


B4  "  FBZR  d"KG  +  -  LC6  -  (•'  FBZR)(GM1} 

tan  BA 

=  (-3414973) (646.8)  -  (24 ,333, 120-3414973) ( 350) 


=  ( -2 . 2088x 109)  -  (7. 32 lx 109 ) 

=  -9.53  x  109 

C4  =  Fbzr  _Lb£  -  LCG  =  3.414973x10^  (172.7) 

=  5.89765837  x  108 

o  =  (1.27857xl07)(3. 414973X106)  -  ( 1, 304378xl06) (5.898xl08) 
(1.304378xl06)(-9.53xl09)  -  ( 1.2785xl07) ( 1.826xl06) 

0  =  4. 4xl013  -  7. 7xl014  =  -7 . 3xl014 

- 1 . 2 10 16  -  2. 3xl013  -1.2x10^ 

=  061  RAD 

3.485 

a-j  =  (cos  .42)  cos  (.262+. 061)  +  .6  sin  (.262+. 061) 

=  (.913)  (.948)  +  .6  (.31646) 

=  .8655  +  1898 
=  1.055324 

by  =  (-cos  .42)  (sin  .323)  +  .6  Lr  (.323) 

=  (-.913)  (.317>  +  .56897 
=  -.2894  +  .56897 
=  +  .279549 

a7 

—  =  1.055324  =  3.775 
/  .279549 

Et  =  3.414973  3.414973 

1  3. 775(cos  .061)  -  (sin  .061)  =  3.775(.998)-.0609 

=  3414973  =  921,203.2 
3.707 

ET 

-  =  921203  =  .946 

T i  io  973700 

•  @  20FPS  Fbzr  =6.829  945xl06  LBS 

A3  =  1  304  378 
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B3  =  1  826  132 

C3  =  Fbrz  =  6.829945  x  106 

A4  =  FBRZ  =  6. 829945xlQ6  =  '2.557xl07 
tanBA  tan  -261 

By  -  Fbrz  d-KG  ♦  LB2  -  LOG  .  (a  .  ^  ^ 

tan  BA 

=  ( -6.829945xl06)  (646.8)  -  (24333120  -  6829945xlQ6)  350 
=  (-4.4176xl09)  -  (6.126xl09) 

=  -1.054371xl010 

C4  =  FBZR  "  LCG  =  6. 829945xlOS(  172. 7) 

=  1,179,531,502 

■■  =  A4C3-A3C4  =  (2.557x107)(6.829x106)  -  (1.3xl06)  1.18xl09 

.  n  n  n"  (1.3xl06)  ( - 1 . 0 54 3x 10 10 )  -  (2. 557xl07)  (1.826xl06) 
H3b4"H4ts3 

=  1 . 7464xl014  -  1 ,534x10* 5 
-1.371xl016  -  4. 669xl013 

=  -1.359xl015  =  .098788  RAD 

-1.3756x1016 
=  5.66 

a?  -  cos  .42  cos  (.262  +  .0988)  +  .6  sin  (.262  x  .0988) 

=  .913  (.936)  +  .6  (.353) 

=  .8542  +  .2118 
=  1.066 

b7  =  (-cos  .42)  (sin  .361)  +  .6  cos  (.361) 

=  (-.913)  (.353)  +  .56136 
=  -.32289  +  .56136  = 

=  .239 

fz.  =  --^-g  =  4.4589 
b7 
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1 


D 


h 


r  _  6.829  x  106 

lT  4.4589Tcos.0988)  -  sin  (.0988) 

=  6.8299xl06  =  6.8299xl06 

4. 4589 (.995)  -  .0986  4.4371-.0986 

=  1. 574xl06 

E._  =  1. 574xl06  =  1.6143 

—  973700 

t1B2o 


POLAR  STAR 

Spread  angle  compliment 


28'  WL 

W.L.  Entrance  =  60°  EA  =  30° 

~2 

Bow  4  =  15° 

B  =  SAC  =  arctan  (sin  15)  =  tan  -1  .25881  =  tan“* 

tan  30  .57735 

B  =  24.146°  =  .42139  Radians  .42 

SA  90  -  24.146°  =  65.8° 

2  SA  =  131.7° 

30'  WL 

W.L.  Entrance  =  59.5° 

Bow  4  =  14.5° 

B  (SAC)  =  arctan  sinBA 
tan  Ea 

T 

B  =  arctan  .2503  =  arctan  .4425 
.5657 

B  =  23.8715° 

B  =  .4165  .42  Radians 

SA  =  90-23.87  =66.13 
2  SA  =  132.26° 


EA  =  29.5° 


=  tan"^  sin!4.5 
tan  29.5 


.44828 
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